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PREFACE 


The data-processing methods and ice data products derived from Seasat radar altimeter 
measurements over the Greenland ice sheet and surrounding sea ice are documented in this first 
volume of a series. The corrections derived and applied to the Seasat radar altimeter data over 
ice are described in detail, including the editing and retracking algorithm to correct for height 
errors caused by lags in the automatic range tracking circuit. The methods for radial adjustment 
of the orbits and estimation of the slope-induced errors are given. The various levels of ice data 
sets are described in this report, but the user is referred to Volumes 2 (Greenland) and 4 
(Antarctica) for more detailed descriptions of the gridded elevation data sets and the 
gco-referenced data bases. 
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SECTION 1.0 
INTRODUCTION 


This volume is the first in a series documenting the data-proccssing methods and ice 
data products derived from satellite radar altimeter measurements over the ice sheets of 
Greenland and Antarctica and surrounding sea ice. The data-processing procedures and 
corrections derived and applied to the Seasat radar altimeter data are described in detail in this 
report. A flowchart depicting the procedures involved in obtaining the various data products is 
given in Figure 1. A detailed description of the editing and retracking algorithm is given in 
Section 2, along with descriptions of the other corrections. The methods for radial adjustment 
of the orbits and estimation of the slope-induced errors are described. The various levels of ice 
data sets produced are described in this report, but the user is referred to Volumes 2 and 4 for 
more detailed descriptions of the gridded elevation data set and the gco-rcfcrenccd data base. 

The input Seasat radar altimeter data, in the form of Geophysical Data Records (GDR's) 
and Sensor Data Records (SDR's) produced by NASA's Seasat project at the Jet Propulsion 
Laboratory, were obtained from the NOAA Environmental Satellite Data and Information Service 
(EDIS) archive on about 1000 magnetic tapes. Development of the data processing methods, the 
production of higher-level geophysical data products, and analysis and evaluation of the data have 
been supported at the Goddard Space Flight Center by funding for research and data analysis, 
provided primarily by NASA’s Ocean Processes Program and by the Climate program. Computci 
programming and technical assistance has been provided by the EG&G Washington Analytical 
Services Center, Inc. until January 1989 and by ST Systems Corporation since then. Numerous 
other individuals have provided valuable assistance. 

Results have been reported in refereed scientific literature (e.g., Brenner ct al., 1983; 
Martin et al., 1983; Zwally et al., 1983; Thomas et al., 1983; and Gundestrup et al., 1986). In 
addition, elevation data in various forms have been provided to other scientists and placed in the 
National Snow and Ice Data Center (NSIDC) and the National Space Science Data Center (NSSDC). 
The purpose of this series of reports is to document technical details and provide guidance to 

users of the ice data products. 

While all reasonable quality-control efforts have been made to eliminate erroneous data, 
some data of questionable quality is likely to have persisted, particularly in the lower-level data 
products. Users should apply normal standards of scientific caution in their use of the data. 
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Figure 1. Processes Involved in Obtaining Data Products 
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The current list of reports is: 


Satellite Radar Altimetry over Ice, Volume 1 : Proces s ing and Corrections of Seasat Data 
over Greenland , July 1989. This volume. 


Satellite Radar Altimetry over Ice. Volume 2 : User's G uide for Greenland Elevation Data 
from Seasat , July 1989. NASA Reference Publication * 

Satellite Radar Altimetry over Ice. Volume 4 : User's G uide for Antarctic Elevation Data 
from Seasat , July 1989. NASA Reference Publication. 

Volume 3 will be the Antarctic equivalent of Volume 1. Additional volumes will include 
descriptions of the data sets being produced by NASA from the radar altimeter data acquired by 
the U.S. Navy’s GEOSAT, using methods similar to those for the Seasat data. 


The Seasat spacecraft (e.g. , Lame and Born, 1982 and Lame et al., 1980) was launched 
in late June 1978, and during its brief 110-day lifetime, collected 90 days of nearly continuous 
radar altimeter data from July 9 through October 10 between the latitudes of 72°S and 72°N. 
Although designed only for measurements over water, the Seasat radar altimeter (MacArthur, 
1978; Tapley et al., 1982; and Townsend, 1980), acquired more than 600,000 useful altimeter 
range measurements over the continental ice sheets of Greenland and Antarctica. 

Over sloping and undulating surfaces, such as ice covered land, or surfaces with 
highly-variable reflecting characteristics, such as in regions of sea ice, the range to the surface and 
the characteristics of the received radar pulse changed faster than the response capability of the 
altimeter electronics. Consequently, it has been necessary to correct each range value for lags of 
the altimeter range servo-tracking circuitry by a procedure called retracking (Martin et al., 1983). 
The retracking correction typically had a mean value of + 1.4 m as applied to the surface 
elevation, a standard deviation of 2.9 m, and maximum and minimum values of ± 15 m. In 
addition, the pulse-limited footprint (1.6 km minimum diameter), which was located near the 
satellite nadir point over the relatively flat ocean, was in general located anywhere within the 
beam-limited footprint (22 km in diameter) over sloping surfaces. The resulting slope-induced 
error, which was nearly 80 m over slopes of 0.8 degree, can be partially corrected using the 
procedures described in Brenner et al., 1983. Corrections are also made for errors in orbit 
determination, atmospheric propagation path-length variations, and earth and ocean tides. 

Elevation measurements were obtained at 0.1-second intervals, corresponding to 662-m 
intervals along the subsatellite ground track. The precision of the corrected range measurements 
is about 1.6 m overall with a minimum of about 0.25 m in the smoothest regions of the ice sheets 
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(Zwally et al., 1983). The 5- to 10- cm precision over the ocean is for 1-sec data averages.) The 
absolute accuracy of the elevations is primarily determined by the limitations on the correction 
methods for the slope-induced errors and uncertainties in the geoid reference level. 

The principal ice data sets produced and/or retained are: 

Level 4: Contour maps and gridded elevations with respect to earth ellipsoid and sea 
level (e.g., this Volume and Volume 2). 

Level 3: Geo-referenced data base including all individual elevation measurements 
(including time, latitude/longitude positions, and slope-correction estimates) accessible 
by geographic cells (e.g., this Volume and Volume 2). 

Level 2: Ice Data Records (IDR's). Orbital-format data records including altimeter 
parameters, corrected elevations, latitude/longitude positions, AGC, applied corrections, 
retracking beta parameters, and estimates of along-track and cross-track slope corrections, 
(this Volume) 

Level 1: Waveform Data Records (WDR's). Orbital-format data records including waveform 
amplitudes by gate, ranges, AGC, and latitude/longitude positions, (this Volume) 

Sensor Data Records (SDR's) 

Geodetic Data Records (GDR’s) 
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SECTION 2.0 
ICE DATA RECORDS 


The Seasat altimeter data were released in two forms: the Altimeter Sensor Data Record 
(hereafter referred to as SDR), and the Geophysical Data Record, GDR. The SDRs were obtained 
from the NOAA/EDIS archives and contain, among other quantities, the telemetered range 
measurements between the spacecraft and earth’s surface, averaged radar return pulses, the 
altimeter status flags and the satellite latitude, longitude, and elevation. The data are output in 
0.098-sec intervals. The GDR’s contain processed SDR data averaged over 1-sec intervals, and 
the sensor, atmospheric, and surface dynamic corrections necessary to utilize the data in detailed 
geodetic work. Data over the ice sheets are not available from the GDR’s. 

To obtain the ice sheet elevation measurements, data from the SDR s are used and the 
appropriate corrections and adjustments applied. This subset of ice sheet data obtained from 
the SDR’s is referred to as ice data records or IDR’s. A detailed description of these records may 
be found in Table 1. The surface heights, located in bytes 73-76 of the IDR, are referenced to the 
IUGG 1980 Geodetic Reference Ellipsoid (Moritz, 1980), which is defined with a 6378. 137-km 
semi-major axis of the earth and a flattening ratio of 1/298.257. Heights relative to sea level can 
be calculated by subtracting the geoid value from the surface height. Geoid values, linearly 
interpolated from a one-by-one degree GEM10-B geoid grid, are located in bytes 61-64 of the IDR. 


Figure 2 is a map of Greenland which depicts the coverage obtained from the IDR’s after 
data were edited and retracked (see Section 2.1). The gaps in the data are a result of the 
altimeter not being able to maintain valid height measurements over the rougher surfaces of the 
ice sheets. Table 2 gives a concise catalog of the available Seasat Greenland IDR data. Included 
in this table are the start and stop locations of each rev, the number of points in each rev, and 
the data base bins (see Section 4.0) through which each rev traverses. The rev numbers are 
ordered such that all ascending passes are listed first, ordered by increasing latitude as they cross 
315 degrees East Longitude. Then the descending passes are listed using the same ordering 
criterion as for the ascending passes. 

2.1 EDITING AND RETRACKING 

As explained in Section 1.0, Seasat altimetry returns over non-ocean surfaces required 
special processing in order to calculate meaningful height measurements. To understand this 
processing one must first have an understanding of the return itself. 
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Figure 2. Seasat Greenland Groundtracks 



Each altimeter return, referred to as a waveform, consists of the output of a set of 63 
gates that span a height window of approximately 30 m. Each gate has a level of return 
associated with it measured in counts. A typical ocean return from Scasat is presented in Figure 
3. The level of return in the first 22 gates is at the noise or pre-pulse level of 4 or 5 counts. The 
level quickly increases to a relative maximum and then slowly decreases over the latter portion 
of the window. There are three half-gates at the center that have a spacing of 23 cm instead ol 
46 cm. The tracking gate is the center of these. The on-board tracker attempts to keep the 
center of the return leading edge positioned at the tracking gate by predicting the travel time of 
each pulse based on previous returns. The measurement telemetered from the altimeter is 
equivalent to the travel time to the tracking gate. 

Altimeter returns over non-ocean surfaces vary greatly from this ocean return. Figure 4 
shows representative returns over ice sheet surfaces for a Scasat pass over Antarctica (Martin 
et al., 1983). The Figure 3 sea ice returns arc represented by one or more sharp spikes that may 
or may not be at the tracking gate. As the altimeter travels onto the ice shelf, acquisition is lost, 
represented as a flat return. On the ice shelf the returns are shaped similar to the oceans, but 
again are not always centered at the tracking gate. As the satellite moves over the ice sheet, 
acquisition is again lost temporarily. Over the ice sheets the returns arc noisy, have multiple 
leading edges, and the mid-point of the first leading edge is not always aligned with the tracking 
gate. 


The measurement telemetered from the on-board tracker needs to be corrected for the 
variation of the mid-point of the leading edge from the tracking gate. This retraeking correction, 
AH ret is calculated as 

AH ret = (Gm-Gt)*g2m (2.1) 

where 

Gm = gate of the mid-point of the leading edge (see Sections 2. 1.3-2. 1.4), 

Gt = the tracking gate (29.5 where the whole gates arc numbered from 0 to 59; sec 

Figure 3), and 

g2m = 


the conversion from gates to meters = .4684375 m/gate. 






It then follows that 


H 


ret. 


H 


+ AH 


ret, 


t+1 


( 2 . 2 ) 


where 


H re t t = the retracked altimeter measurement at time t, 

H mea s t — the measurement calculated by the on-board tracker at time t, and 

AH fCt ^ = the retracking correction calculated from waveform at time t + .098 see. 

Due to the return being telemetered one time step later, the retracking correction for the 
measurement at time t is calculated from the return at time t + .098 sec. Methods have been 
developed at NASA/GSFC to calculate the AH ret for returns over the ice sheet, ice shelf, and sea 
ice which can yield valid height measurements. A detailed description of these procedures may 
be found in Sections 2. 1.3 and 2. 1.4. Parameters resulting from these retracking techniques may 
be found in bytes 109-144 of the IDR. The criteria used to automatically select and discriminate 
between different types of returns are described in the next two sections. 

2.1.1 Selecting Retrackable Non-Ocean Altimetry Returns 

The SDR for Scasat includes all telemetered altimeter data even when the instrument 
was in calibration or standby mode. Since valid measurements could be acquired when the 
tracker was in acquisition mode, all data that are not in acquisition or track modes are discarded. 

All tracking and acquisition returns have to meet two initial tests to determine if the 
waveform actually represents the initial return, or if the return is outside the tracking window. 

1) The counts in the first gate must be less than 100; 

2) There must be at least one gate with a count value greater than 25. 

2.1.2 Categorizing the Returns 

The remaining returns are then categorized into two groups. Group one will be referred 
to as specular and consists of those returns that display a sharp spike. Returns in this category 
are usually found in regions of sea ice or over flat, desert-type surfaces. The second group, 
consisting of the remaining returns, is called diffuse and resembles ocean returns. These returns 
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are found over continental ice and the ice shelves. Different methods are used to retrack each 
group. 


Returns are automatically categorized as either diffuse or specular depending on the 
existence of a significant spike in the return. To determine this the following algorithm is used. 
The noise level, Yn, is calculated as the average number of counts in the first five gates. The 
maximum, Ymax, is calculated as the maximum number of counts in any gate. The value Ymed 
is then calculated using the equation 

y mcd , IVtmstM + Yn (2.3) 

2.0 

The gate number, Gmid, is then found as the first gate where the number of counts exceeds 
Ymed. Two sums of consecutive counts from the signal are then formed, Ylow and Yhigh, where 


i=Gmid + 9 
Ylow 2 Y, 
i=Gmid 


(2.4) 


Yhigh 


1 = Gmid+20 

2 Y, 

i=Gmid + 10 


(2.5) 


If Gmid is so large that there are less than 20 remaining gates, then the number of gates used 
to form the sums is adjusted. When the ratio of Yhigh/Ylow is < 0.7, the return is considered 
specular. 

2.1.3 Retracking Specular Type Returns 


Specular waveforms are not found in the Seasat altimeter data over Greenland. This is 
probably due to the absence of sea ice near Greenland during Seasat s lifetime. As a result, all 
of the Greenland returns are retracked using the diffuse method. However, for the sake of 
completeness, the method used to retrack specularly shaped returns, which is employed in the 
region of the Antarctic, will be discussed. 

Specular-type returns are defined for this procedure as being characterized by one or 
more extremely sharp spikes and are retracked by attempting to locate the mid-point or half- 
power point of the first significant spike. In addition, since the shape of the return essentially 
records topographic characteristics, parameters are also calculated which define the shape of a 
single-or double-peak return. Figure 5a shows the five-parameters required to define a single- 
peak return, while Figure 5b shows the nine-parameters required for a double-peak return. 
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2.1.3. 1 Half-power Point of First Significant Peak 


In determining the mid-point of the first significant spike, the location of this spike must 
first be found. The value of Ymed, which is calulated to determine whether or not the return is 
specularly shaped (Equation 2.3), is used. Starting with the gate number prior to Gmid, where 
Gmid is define to be the gate number whose counts exceed Ymed, a gate is sought whose counts 
exceed or equal 25% of the difference between Ymax and Yn. Upon finding this gate, Grise, it is 
determined to be the first significant spike if the following conditions are met: 


Y 


Grise + 1 


Y 

1 Grise 


< 0 . 


for Y Grlse > Ymax * 3 


(2.6) 


where 


Y Crlse is the counts for gate Grise, and 

^Grise+l * s counts f° r gate Grise -I- 1 . 

Smaller, more rounded waveforms, which might be encountered in the vicinity of an ice 
shelf require that the following condition be met: 


Y 


Grise + 1 


1 Grise 


< (Ymax-Yn) * .05 


for Y Grlse ^ Ymax * 3 • 


(2.7) 


Grise is incremented by one, up to the maximum number of gates, until one of the above 
conditions is met, after which the gate of the first significant spike, Gist, and its corresponding 
counts, Ylst, are used to determine the half-power point of the peak. The count value at the half- 
power point, Ymidl, is determined as follows: 

v (Ylst-Yn) , „ 

Ymidl = + Yn . (2.8) 

2.0 
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The exact gate location of the half-power point, Gtmidl, is then determined by performing 
a linear interpolation for the count value Ymidl located between gates XI and X2, with 
corresponding count value Yl, Y2. 


2. 1.3.2 Remaining Parameters to Define Shape 

In order to define the exact shape of the specular returns depicted in Figures 5a and 5b, 
it is necessary to calculate several other parameters in addition to the noise level, the maximum 
counts of the first significant peak, and the gate location of the half-power point. For the single- 
and double-peak return, additional quantities which define the width of the significant peak and 
slope at the half-power point are defined. A double-peak return has four additional quantities 
calculated: the maximum counts for the second significant peak, the gate location of the half- 
power point for the second peak, the slope at the second half-power point, and the minimum 
counts found between the two significant peaks. 


The slopes at the half-power point for both the first and second significant peaks. Sip 1st 
and Slp2nd, are determined by the following algorithm: 


Sip 1st 


Y2 - Yl 
X2-X1 . 


(2.9) 


Slp2nd uses the gate locations and corresponding counts determined to surround the half-power 
point of the second significant peak. These values are found in a manner similar to that of the 
first peak. 


The actual existence of a second significant peak is determined in the following manner. 
Starting with the gate location of the first significant peak, the difference between counts of 
consecutive gates is monitored. As soon as the change in successive gates becomes negative, at 
gate location Gcntmin, it is assumed that another peak has been encountered. At this point, a 
sum if formed, Totup, which totals the counts in all gates following the Gcntmin. When Totup 
equals or exceeds 9% of Y 1st then the second peak is considered significant. The gate at which 
the second peak occurs, X2nd, is determined to occur when the difference in the counts of 
consecutive gates becomes positive. 


The counts at the second significant peak, Y2nd, are then used In the following manner 
to calculate the counts at the half-power point of the second peak, Ymid2: 

Ymid2 = (Y2nd-Cntmin) + Cntmin (2.10) 

2 
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Again, a linear interpolation is performed in a manner identical with the first significant 
peak to determine the exact gate location of the second significant peak half-power point, Gtmid2. 

The final parameter to be determined is the total width of the peak or peaks at the first 
hair-power point. The width is defined as the number of gates between Gtmidl (Section 2. 1.3.1) 
and the location, Gtrail, where the trailing edge passes through Ymidl (Equation 28). The width 
is computed as follows: 

Width = Gtrail - Gtmidl . (2.11) 

In summary, the parameters for a specular return with a single significant peak are as 

follows: 


01 = Yn 

p 2 = Ylst 

(3 3 = Gtmidl (2.12) 

0 4 = Slplst 

P 5 = Width . 

The parameters for a specular return with double significant peaks are as follows: 

01 = Yn 

)8 2 = Ylst 

/3 3 = Gtmidl 

04 = Slplst 

fi 5 = Y2nd (2.13) 

= Gtmid2 
P 7 = Slp2nd 

P 8 = Width 

(3g = Cntmin . 

2.1.4 Retracking Diffuse-Type Returns 

The method used to retrack the diffuse return is to model the return with a function that 
has the retracking position (the mid-point of the leading edge) as a parameter. The Bayesian least- 
squares method (Ref. 8) is used to solve for the parameters of the function that best fit the 
return. For this method, initial estimates of the parameters must be provided. Weights are given 
to these initial estimates that designate how well each parameter is known relative to the others. 
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Residuals are then calculated between the return value and the function value at each gate. These 
residuals are weighted based on their proximity to the mid-point of the leading edge position. A 
minimum to the sum of these squared weighted residuals is sought by an iterative method which 
simultaneously adjusts all of the function parameters. The process is repeated until convergence 
or until the maximum number of iterations is reached. Because linear methods are used to solve 
a non-linear problem the procedure can be numerically unstable. Checks are done to assure the 
reasonableness of the results. The key to making this method function correctly is in the choice 
of the initial estimates and weighting functions. 

The theory of solving for the function parameters using Bayesian least-squares can be 
found in Ref. 8. The actual equations used will be presented here without justification. 

Given an overdetermined set of equations MX=R where 


M = the matrix of partials 


dc t dc, 

#1 ' ' ' 




m > n 


( 2 . 14 ) 


x = column vector = 


fid - fii 


Pen ~ Pn 


( 2 . 15 ) 


R 


m i -°l 


m 


m 


( 2 . 16 ) 


and 

m t = observed value (counts at t=gate i) t 

c L = calculated values of nij based upon a given set of parameters /?, 
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j3j = current best estimate of the model parameters /3, 
£ c j = corrected best estimate of the model parameters /3, 
i = gate number (0 - 59), and 
n — number of parameters in the function. 

We can then define a weight matrix, W 


wtj 0 


W = 


0 wt 


(2.17) 


where wt { is the weight associated with each observation i. 

If we multiply both sides of the equation by W we get 


WMX = WR . 


Multiplying through by M T gives 

M T WMX=M T WR . (2.18) 

The solution of X is solved for as 

X = (M T WMr 1 M T WR (2.19) 

where M T WM is referred to as the normal matrix. To add information as to the validity of the 
current best estimate of the model parameters the a priori covariance matrix V Q is included 


Vo = 


wt/31 0 


0 wt 


-/3n 


( 2 . 20 ) 


17 



where wt^j = weight associated with the a priori value of parameter j. This matrix is then added 
to the normal matrix before it is inverted so the equation becomes 

X = [M T WM+V 0 ] -1 M T WR . (2.21) 

X then is the vector giving the new best estimate of the /3 parameters. 

2. 1.4.1 The Function Representing the Altimeter Return 

It has been shown (Miller and Brown, 1974) that the mean return waveform over a 
Gaussian surface can be mathematically described using the function 

c(t) = jSj 4 P 2 * P (W) (2.22) 


where 

P(W) 


W 

J Z(q) dq 

- X 




(2.23) 

(2.24) 

(2.25) 


This assumes that the pointing angle errors have negligible effects on the waveform shape. This 
also represents the ice sheet waveforms very well if it is modified to include a slope to the trailing 
edge. The modified function used to represent the diffuse-type waveforms is chosen as 

c(t) - /3 1 +j3 2 (l+j3 5 Q(x)) P (W) (2.26) 

where 

Q(x) = 0 for t < + 0.5 

= t-x for t > H- 0.5 P 4 
This is plotted in Figure 6a where 
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04 and 37 are risetime parameters for the 1st and 2nd ramp respectively 
(t -83) (t - 36) 

Where y = 3i + 32'P (1 + 3gQ( x i)) + (3sP — ~ — (1 + 33(Q(^2» 

34 37 


Xi = t- 33 -0.5 3 4 Q(x) = 0 for x < 0 p( Z ) _ 

X 2 = t - 3 6 - 0.5 3 7 = 1 for x _> 0 

Figures 6a and 6b. Diffusely Shaped Waveforms 


P(z) = ( —=• exp (- q 2 /2) dq 
_ « v2n 





x = 03 + 0 5 ^4 ’ 


The partials of this function with respect to each parameter are 

*L = l.o 

= P[W] + /3 5 9 P[W] 

^2 

= -£ 2 {- n+/3s9) ^L+ P(W) /3 5 } 

<P 3 ^4 aw 

*£_ = / 3 2 { (1+/3s9) «Lw^ 5 PtW) } 

^4 /3 4 dW 2 

= /3 2 Q P[W| 

^5 

where 

!L = _L exp (iW 2 ) . 

aw V2 tt 2 


( 2 . 27 ) 


( 2 . 28 ) 


( 2 . 29 ) 


( 2 . 30 ) 


( 2 . 31 ) 


The value of /3 3 is the mid-point of the leading edge, Gm. As previously noted, some of the 
returns display multiple leading edges. A nine-parameter function is used to represent these 
returns, where the mid-point of the first leading edge is still £ 3 . The mid-point of the second 
leading edge, f36, probably represents a return from another surface in the footprint and is being 
stored for future use. The nine-parameter function is 

c(t) = fl 1 +/3 2 P(Wj) (1 +/3 9 9(x 1 )) + fi 5 P(W) (l+/3 8 (Qfx^) (2.32) 

This is plotted in Figure 6b where 

X! = H9 3 -5 £4 

X2 = t-^g-0.5 J3 ? 
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04 

w 2 = 

P 7 

The partials of this nine-parameter function are 


II 

1.0 

(2.33) 

51* 

11 

P(W,) [l+jSglQ! 

(2.34) 

II 

-P 2 [ (1+/399l) £ + P(W j) fig] 
fi, «W, 

(2.35) 

51* 

II 

_o piw.ift, , u^ 9 e L ) 3P Wi , 

2 p A aw l 

(2.36) 

II 

1 + 08 9 2 P(W 2 ) 

(2.37) 

ac_ = 

-/3 S [P(W 2 )/3 8 +lil0£ 9 2 ) *L] 
0 7 dw 2 

(2.38) 
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dc = 

— = @5 Q 2 p ( W 2> (2.40) 

= P 2 92 P ( W l) • (2.41) 

$9 

2. 1.4.2 Setting the Initial Estimates for the Parameters 

Initial estimates of each parameter are calculated from each individual return. To 
calculate these the general shape of the waveform is mathematically described by defining a mean 
slope and average value (bias) for every whole gate. For gates 4 through 56, the mean slopes and 
biases correspond to a straight line that is fit using least-squares minimization through the gate 
in question and the six surrounding gates. The biases for gates 1 through 4 are taken as the 
gate values and the slopes arc defined as zero. For gates 57 through 60 the biases are the gate 
values and the slopes are defined as the slope calculated for gate 56. This set of slopes and 
biases is then interrogated to determine the locations of the leading edges and how many occur 
in the waveform. 

The conditions required for a leading edge at gate Ir are: 

1) The Slopc(Ir) must be greater than a given value, Thsl. A value of Thsl = 0.5 count/gate 
is used to find the first leading edge, for succeeding leading edges Thsl is set to 1.0 
count/gate. These numbers were chosen by visually and mathematically evaluating 
many typical ice sheet waveforms to determine when a leading edge designating a valid 
return could be perceived. 

2) The Slope(Ir) must be a relative maximum, ic: 

Slope(Ir) > Slope (Ir— 1 ) 

Slopc(Ir) > Slope (Ir+1) . 
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3) There must be a significant increase in counts after the leading edge compared with 
that before the leading edge, i.c.: 

Bias(Ir4-3)-Bias(Ir-3) > Thbs 
where 

Thbs = 13.5 counts for first leading edge 

= 20.0 counts for succeeding leading edge. 


4) If there was a leading edge already detected within 3 gates of Ir then the location is 
taken as that with the larger slope. 

The initial estimates of the function parameters are then calculated from the position of 
the leading edge(s) and the Slopes and Biases. The five-parameter function (2.26) is used when 
only one leading edge is found, the nine-parameter function (2.32) is used when two or more 
leading edges are found. 


Initial estimates, and the corresponding standard deviations of these estimates, Sig(l) 
through Sig(5), for the five-parameter function arc defined as: 


0 ? 

02 

03 

04 


05 


Bias(4) (counts) 

Sigd) 

= 0.01 (count) 

Bias(Ir + 3)-Bias(4) (counts) 

Sig(2) 

= 10.0 (counts) 

Ir (gate) 

Sig(3) 

= 1 0 O (4) (gates) 

(|Bias(Ir + 3)-Bias(Ir-3)l/ 

Sig(4) 

= .01 /3 {) (4 ) (gates) 

Slope(Ir)}*0.5 (gate) 
0.0 (count/gate) 

Sig(S) 

= .01 (count/gatc). 


(2.42) 


Initial estimates and the corresponding standard deviations for the nine-parameter function arc 
defined as: 


p° 

= Bias(4) (counts) 

Sig(l) 

= .01 (count) 

P°2 

= Bias(Ir 1 43)-Bias(4) (counts) 

Sig(2) 

= 0.1 (count) 

PI 

= Irl (gates) 

Sig(3) 

= .05 /3 0 ( 4) (gates) 

P°4 

= {[Bias(Ir 1 +3)-Bias(Ir 1-3)1/ 
Slope(Irl)}*0.5 (gates) 

Sig(4) 

= .005 )3 0 (4) (gates) 

PI 

- Bias(Ir2 + 3)-Bias(Irl +3) (counts) 

Sig(5) 

= 0.1 (count) 

PI 

= Ir2 (gates) 

Sig(6) 

= ,05/3 o (7) (gates) 

P°7 

= {[Bias(Ir2 + 3)-Bias(Ir2-3)] 
Slope(Ir2)} (gates) 

Sig(7) 

= .005 /3 q ( 7) (gates) 

'CO 

oco 

= 0.0 (count/gate) 

Sig(8) 

= .01 (count/gate) 

p°9 

= 0.0 (count/gate) 

Sig(9) 

= .01 (count/gate) 


(2.43) 

where 

Irl is the predicted gate corresponding to the mid-point of the first leading edge 
Ir2 is the predicted gate corresponding to the mid-point of the second leading edge. 

2. 1.4.3 Calculating the Weight Matrix, W 

The weight associated with each observation, wt t , is selected to optimize the fit in the vicinity 
of the leading edge. 

wt, = 1 + Kj ♦ [exp (K 2 ) + K 3 ] (2.44) 

where 

K, = (Iter" D * 0.5 

I ter = iteration number 

K 2 = T c + 0.5 

= Min (K 2 , 60) 

= Max (K 2 , 1) 

T c = X { - /3 3 - Max (5.0, /3 4 ) for 5-parameter function 
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= X[ -y3 g - Max (5.0, /3 7 ) for 9-parameter function 
X, = gate number of the ith observation 

for the five-parameter function 
K 3 = 0 for |T C | 5: 2.0 

= 1 for |T C | < 2.0 

for the nine-parameter function 
K 3 = 0 for |T C | ;> 5.0 

= 1 for |T C | < 5.0 . 

2. 1.4.4 Calculating the Covariance Matrix, V 0 

A priori values of V Q are calculated from the sigmas in equations (2.42) and (2.43) as 

follows: 

W t = mscale/Sig(j) 2 (2A5) 

P} 

u> scale = 1 + .6 * K * Hl/3/( 120*g2m) (2 46) 

Hl/3 = 1.875 * /3 4 

K =4 

Using the function, m scale, causes the initial estimate information to have a greater effect 
on the solution when the rise time is large. 

After each iteration, n, the values of Sig(3), Sig(4) and K are altered as follows: 

Sig(3) = SigOl^j * 0.1 

Sig(4) = Sig(4) n _! * 10.0 
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+ .5 . 


K = K n-l 

This has the effect of weighting the current best estimate of the leading edge position 
more and the rise time of the leading edge less. This has proven to speed up convergence. 

2. 1.4.5 Method of Iteration 

An interative scheme is used starting out with the initial estimate of the (3 parameters. 
The Bayesian least-squares method is then used to solve for another set of /3 parameters that 
better fits the data. Iterations are performed always using the current set for the best estimate 
until AH rct , as calculated from /? 3 (2.1), converges to within 10 cm or the number of iterations 
exceeds 7. 


Each succeeding set of (3 parameters is checked for reasonableness using these criteria: 

0.0 < p 2 

0.0 < p 3 < 60.0 

0.0 < 0 4 

(3 3 < < 60.0 

o.o < p 7 . 

If any of the criteria fail, then the fit is considered unsuccessful and the waveform is discarded. 


After convergence or the maximum number of iterations is reached, tests are then made 
to assure that the values reasonably represent the return. The rms of the residuals between the 
waveform and the function for the portion of the waveform from gate zero to just past the top of 
the leading edge is calculated. 


Iedit 

2 (C t -m ( ) 2 

RMS e = A= l 

Iedit 

where 

Iedit = |S 3 + 0.5 /3 4 for the five-parameter function 
= f3 6 -I- 0.5 (3 ? for the nine-parameter function. 


If RMS e is greather than 20.0 counts then the fit is unacceptable. If the nine-parameter 
function is being fit and the process is unsuccessful, then the initial estimates are reset to 
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coincide with the initial estimates for the first leading edge and a five-parameter fit is tried. II 
problems occur during the five-parameter fit, the initial estimates are altered so that the leading 
edge position is taken as the gate, I r , where Slope (I r ) (as defined in Section 2. 1.4.2) is a 
maximum for the waveform. If the fit is still unsuccessful, then the waveform is discarded. 

The procedures explained here and the numerical values given yield the best results to 
date. Wherever possible values were chosen based on theory, but many times trial and error was 
necessary. At the time the Seasat Greenland data were processed, the procedures and numerical 
values differed slightly. There was no RMS E check as explained in the last part of Section 2. 1.4.5, 
nor were the initial parameter values altered if an unsuccessful fit was made. The variables that 
were different and their values for the Greenland processing were: 

Thbs = 5.0 counts for the first leading edge 

= 10.0 counts for the second leading edge 

Sig(3) = /3 Q (4) (for the five-parameter function) 

Sig(4) = O.l/3 0 (4). 

A direct consequence of these differences was that the entire Greenland data set had to be visually 
reviewed to assure that the fit adequately represented the data. This resulted in approximately 
1% of the data being discarded which would not have been rejected using newer methods. The 
newer methods described here identify these problems automatically. 

2.2 SENSOR-RELATED CORRECTIONS 


After the ice altimeter data are edited and retracked, the precise orbits from NASA/GSFC 
(PGS-S4) are used to calculate the measured ice sheet elevation above the ellipsoid (Lerch et al., 
1982). Corrections are then applied to correct for sensor-related biases. 

Both the time tag and center of gravity corrections are calculated using the algorithms 
released by JPL (Lorell, 1979). These are summarized below. 

2.2.1 Time Tag Correction 

The SDR time tag, t SDR , is corrected for a track mode correction and a signal travel time 
correction so that the resultant data time, t, refers to the time of signal reflection from the ice 
sheet. 
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- 0.0794 + H/c 


(2.47) 


4 - ^SDR 

where 

c = 2.99792458 x 10 8 m/sec, 

H = spacecraft altitude in meters, and 

0.0794 is the track mode correction in seconds. 

2.2.2 Center of Gravity Correction 

The correction applied to make the spacecraft center of gravity the height reference point 
is 


AH cg = 


Z cg" 


(2.48) 


where 


Z cg ~ the distance from the altimeter base plate to the spacecraft center of gravity. 

This varied during the flight due to maneuvers. Table S-07 of Lorell (1979) 
is used to obtain Z fHr 

^cone — “ 1.238 m which is the sum of the distance from the feed flange on the antenna 

to the base plate and a distance corresponding to a time bias in the electronic 
circuitry. 

This correction is located in bytes 49-52 of the 1DR. 

2.3 ATMOSPHERIC CORRECTIONS 


The measurements are corrected for ionospheric and tropospheric refraction using 
parameters supplied by JPL on the GDR’s (Lorell et al., 1980). 

2.3. 1 Ionosphere Correction 

The ionosphere correction for the ice data, AHjq^, is calculated by linearly interpolating 
from the ionosphere corrections on the GDR's. Bytes 57-60 on the IDR contain the value of this 
correction. A detailed description of the algorithm used is given in Lorell et al., (1980). 
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2.3.2 Troposphere Correction 


The wet tropospheric correction is calculated using the following equations explained in 
Lorell et al., (1980). 


AH TROPwet = 2.277*10-^(1.25503^+0.5) 


(2.49) 


where 


(7.5*T r -273. 16)/(T k -35 .86) 
E q = 6. 1 1*H r * 10 K * 


T k is the surface temperature calculated by assuming a linear temperature profile with 
boundary conditions: 


at sea level T K = 2. 73. OK 


at 3200m above sea level T K = 243. OK, and 


H r is the relative humidity (assumed to be 100% over the ice sheet). 


The dry tropospheric correction is calculated from the equation 


AH. 


TROPi 


DRY 


= 2.277* 10 _3 *{P* [ 1.0+0.0026*cos(t£)]} 


(2.50) 


where 


4 > 

p 

P o 

Ht 


= subsatellite latitude, 


= P Q * (1.0-1.1138*10 


-4 ♦ 


Ht), 


- is the atmospheric pressure interpolated from the GDR’s, and 
is the ice sheet elevation above sea level in meters. 


The total height correction due to the troposphere is 


^TROP - AHtrop,!^ + AH TROP 


WET 


DRY' 


(2.51) 


The troposphere correction may be found in bytes 53-56 of the IDR. 
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2.4 SURFACE DYNAMIC CORRECTIONS 


The solid earth tides are computed by linearly interpolating their values from the GDR’s. 
The resultant interpolated value may be found in bytes 83-84 of the IDR. 

2.5. ORBITAL CORRECTIONS 


The NASA/GSFC PGS-S4 orbits which are used to improve the height measurements, have 
rms radial errors of 1.5 m. In an effort to reduce the radial error of these orbits, a technique was 
devised to further improve the orbit accuracy by referencing the orbits to a common ocean 
surface. Previous attempts to adjust the orbits using crossover minimization techniques with the 
ice sheet crossovers proved unsuccessful due to extreme segmentation of the data (see Figure 2). 
The new technique is not dependent upon the ice data but upon ocean altimetry, and utilizes the 
smoothed Seasat 84306 global ocean surface (Marsh et al., 1986). Through crossover 
minimization techniques the radial orbit error for the 84306 ocean surface has been reduced to 
1 1 cm in the open ocean areas. 

The method involves obtaining the residuals between the Seasat ocean data for passes 
which traverse Greenland, and the smoothed 84306 ocean surface. Using least-squares 
minimization, these residuals are then fit to a linear or quadratic function depending on the 
proximity of the data to Greenland. The function is, in turn, interpolated or extrapolated to 
determine the value of the orbit adjustment over Greenland which is to be subtracted from the 
surface height. This function is of the following form: 

l*(t) = Cq + C j At + C 2 At 2 ^2 52 ) 


where 


C 0 'Ci.C 2 are the coefficients of the fit where the units are meters, meters/fractions of 
a day and meters/( fractions of a day) 2 , respectively, and 

^ is the time from the start of the pass in fractions of a day. 

Since this method attempts to adjust for orbit error only, the ocean data which are used must 
have all sensor, atmospheric, and surface dynamic corrections applied. The ocean data used in 
the adjustment are obtained from the Seasat Geophysical Data Records (GDR’s), as corrected by 
JPL (Lorell et ah, 1980). 
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Since the orbit error is strongly periodic, with a dominant frequency of two cycles per 
one revolution, only data from the northern hemisphere need to be used in computing the orbit 
adjustment over Greenland. 

The distribution of the data aflects the way in which the residuals are fit. To aid in 
categorizing the distributions of data, the northern hemisphere is subdivided into five ocean 
regions: 1) the area to the east of Greenland and within 1000 km. of the coast; 2) the area to 

the east of Greenland from 1000 km. from the coast to the Greenwich meridian; 3) the Indian 
Ocean; 4) the area to the west of Greenland between Greenland and North America; and 5) the 
Pacific Ocean (see Figure 7). The type of fit performed depends upon particular regions 
containing a minimum amount of data. If the criteria are not met, then no fit is performed. 

Figure 7 summarizes the type of fit which is performed depending upon the region(s) in 
which data are found. An ’X’ in regions 1, 2, 3 or 5 represents a minimum of 10 points, while 
region 4, due to its limited open ocean area, requires a minimum of 19 points. Linear fits are 
performed when data are found either very close to Greenland or are widely separated from 
Greenland. Quadratic fits are performed when the data are more evenly distributed over several 
regions. 

After the coefficients for the fit are initially determined, outlying data which satisfy the 
following criterion are removed: 

[H(t)-rtt) | ;> m * RMS (2.53) 


where 


m is an integer editing multiplier, 

RMS is the rms between the residual heights and the function f(t), and 
H(t) is the surface elevation of the datum point. 
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Figure 7. Orbit Adjustment Regions and Effects of Data 
Distribution on the Orbit Adjustment Fit 
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The remaining data are then used to solve for the function. This process is repeated until 
either the latest computed rms does not change by more than .02 m from the previous iteration, 
or 15 iterations are completed. In the case of the Seasat Greenland data, an editing multiplier 
of 4.0 is used with an initial rms of 20.0 m. 

After solving for the coefficients and removing outliers, the function must satisfy a final 
test. For a linear function, the orbit adjustments are computed at the endpoints of the pass. If 
the absolute value of the orbit adjustment at either endpoint exceeds 3.0 m, then the function 
is not used. In the case of a quadratic function, the extremum of the function is first located. 
If the extremum is outside the endpoints of the data just fit, then the endpoints of the pass arc 
checked as in the linear case. If the extremum lies between the endpoints, its value is checked. 
Again, a 3.0 m adjustment is deemed too large and if exceeded, an attempt is made to refit the 
data with a linear function. Of the 331 GDR passes for which an orbit adjustment was computed, 
181 resulted in quadratic fits and 150 in linear fits. Of the 194 quadratic fits initially attempted, 
12 failed the extremum test and were refit using a linear function. Of these, only one failed the 

endpoint test. 

Two examples of results from the orbit adjustment procedure are shown in Figures 8 
and 9. In the first case (Figure 8), data which are found in close proximity to Greenland are fit 
by a linear function. The latitude and east longitude of the points along the pass closest to the 
west and east coasts of Greenland are indicated. A linear function is fit to the smoothed ocean 
surface residuals. The orbit adjustment in the region traversing Greenland is indicated by dashes. 
Figure 9 shows the orbit adjustment results when a quadratic fit is necessary due to data being 
available just off Greenland’s east coast and in the Pacific Ocean. The final rms between the 
data and function are 27 cm in the linear case and 14 cm the quadratic case. 

Table 3 summarizes the orbit adjustments computed for each GDR rev at 310, 320, and 
330 East Longitudes, representing the west coast, central region, and east coast of Greenland. 
Also included are the coefficients for the function (Equation 2.52) and the elapsed time in 
fractions of a day from the start point of the pass used to compute the adjustment for the 

longitude in question. 

Utilizing Equation (2.52), the orbit adjustment is then computed for each Seasat IDR, and 
subtracted from the surface height. The orbit adjustment and its corresponding rms are located 
in bytes 93-96 and 97-100, respectively, of the IDR. 
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Application of the orbit adjustment to the data yields improved crossover results. When 
the differences in heights are computed at 1235 crossover locations for ascending and descending 
passes over Greenland, the resultant crossover residual mean of the data without the orbit 
adjustment is 33 cm with an rms of 1.15 m. After application of the orbit adjustment, the data 
give a crossover residual mean of 7 cm and an rms of 0.99 m. 

2.6 SLOPE CORRECTION 

The altimeter height is measured to the closest point within its footprint, which does not 
correspond to the subsatellite location for sloping surfaces. This effect introduces an error into 
the height measurement which can be corrected by adjusting either the value of the measurement 
or its location (Brenner et al., 1983). Upon examination of both techniques, the method which 
was chosen for the Seasat data is to adjust the measurement. The magnitude of the slope- 
induced error may be represented by: 

^SLOPE — H( 1-cos a) (2.54) 


where 


H is the satellite altitude in meters 
a is the maximum regional surface slope in radians 


or 

^SLOPE = * for small a. (2.55) 

2 

The surface slope in Equation (2.55) for any one point is calculated using the following 
equation: 


a 



2 

along-track 


+ 


2 

a 

cross-track 


(2.56) 
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where 


a along-track 


is the slope of the surface in the along-track direction of the data, and 


is the slope of the surface in the cross-track direction of the data, 
perpendicular to the along-track direction. 


The cross-track slope is obtained by using a reference surface of Greenland, generated 
from the Seasat data. This surface consists of a two-dimensional grid of heights. The spacing 
between grid points is 20 km. Bilinear interpolation between these grid values is used to 
determine the heights at the points where the cross-track intersects the closest grid lines. From 
these heights, the cross-track slope is then determined. 

The along-track slope is obtained using the available along-track data. Since the height 
profile is initially unknown, an iterative procedure is used to attempt a reconstruction of the 
true height profile. The initial along-track slope at a data point location is calculated by 
performing a linear fit to the live elevations of the along-track data points nearest the data point 
in question. A slope correction is then calculated for that point and each point in the pass using 
Equation (2.55), but applying only 25% of the correction to the elevations. This entire procedure 
is repeated using the revised elevations three more times, each time applying 25% of the current 
elevation correction. After the final iteration, the total along-track height correction and Equation 
(2.55) are used to calculate an "effective" along-track slope. This slope may then be used in 
Equation (2.56) along with the cross-track slope to calculate the total slope. In the case of both 
the along and cross-track slopes, a maximum of .8 degree is allowed. This is a limitation set by 
the physical characteristics of the altimeter. 

If two points cannot be found on both sides of the point being adjusted, after having 
searched 10 km in both directions, then the reference grid which is used to calculate the cross- 
track slope is also used to determine the along-track slope in a manner equivalent to the cross- 
track slope calculation described above. 


Slope corrections are not applied to the surface heights on the IDR's. However, the along- 
track and cross-track slopes, from which the slope correction may be computed, are stored in 
bytes 85-86 and 87-88, respectively. Bytes 89-90 contain the size of the window required to find 
the five points to perform the along-track linear fit. Bytes 91-92 give information pertaining to 
how the along-track and cross-track slopes were determined. 
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2.7 SUMMARY OF CORRECTIONS 

In order to obtain a corrected surface elevation relative to sea level with the solid tide 
effects removed, the following algorithm is used. 

h cor = H sc ~ h alt “ ah ret " Ah cg + Ah ion + Ah trop ” Ah tide 

“ AH ORB “ AH SLOPE “ H GEOID (2.57) 

where 

Hgc is the height of the spacecraft above the ellipsoid, 

H ALT is the ori ginal altimeter measurement, 

A ^ret i s rotracking correction, 

AH cg is the center of gravity correction, 

A ^ION * s the ionospheric correction, 

a H T rop is the tropospheric correction, 

Ah tide is the value of solid tide, 

Ah orb is the orbit adjustment, 

a H SL ope * s the siop 6 correction, and 

h GEOID is the value of the geoid. 

The surface elevation on the IDR is relative to the ellipsoid and is corrected for 
tropospheric and ionospheric effects, the center of gravity offset, the retracking correction, and 
the orbit adjustment when available. However, the elevation still contains solid tide effects, and 
the application of the slope correction or removal of the solid tides have been left to the discretion 
of the user. The surface elevation status word located in bytes 77-78 of the IDR should be 
checked to verify whether or not corrections have been applied. 
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Corrections which are applied to the altimeter measurement are done in the opposite 
sense from the surface elevation corrections and may be verified using the altimeter measurement 
status word in bytes 13-16 of the IDR. 

An outline of the adjustments and corrections required to the Seasat data and their values 
or range of values is given in Table 4. 
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SECTION 3.0 

WAVEFORM DATA RECORDS 


The averaged radar return pulses contained in the SDR’s are stored on a separate file 
called the Waveform Data Records (WDRs) to facilitate their use. Table 5 outlines in detail the 
format of this record. 

The time, geographical position, and altimeter measurement on the WDRs are not identical 
to the corresponding records on the IDRs. This is due to the fact that the WDRs information is 
obtained directly from the SDR's without the application of any correction or adjustment of any 
kind. The time differs by the time tag correction described in Section 2.2.1. Positions on the 
WDR are from the orbits on the SDR’s and not PGS-S4 orbits. The altimeter measurement 
represents the raw observation on the SDR without any of the corrections described in Section 
2.7 applied. 
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SECTION 4.0 

GEO-REFERENCED DATA BASE 


Ordering the Seasat data merely by time presents certain limitations when only data in 
a particular locale are desired. This situation arises when data are used to generate a grid of 
smoothed surface heights. To circumvent this problem, a data base was developed which orders 
the Seasat data by geographical areas or "bins". Figure 10 shows the configuration of the 4,300 
bins in the vicinity of Greenland. Bin sizes vary in order to compensate for the higher data 
density near Seasafs maximum latitude. Each bin is assigned a number starting with "1" in the 
southwestern-most corner. Bin numbers increment first from west to east and then from south 
to north. The ending bin number for each row is indicated in the right-most margin of the map 
in Figure 10, while the number of data points is printed within the appropriate bin. Bins which 
contain no data have no number entered. Table 6 summarizes the number of points and the rev 
numbers found in each bin, along with the geographical coordinates of the southwestern-most 
corner of the bin. The bin number in which a particular data point is located may be found in 
bytes 153-156 of the IDR. 

The geo-referenced data base is a subset of the IDR's, containing only information relating 
to the position, rev number, surface height, slope correction and orbit adjustment for each data 
point. Slope correction and orbit adjustment values are flagged with a -9999, if unavailable. In 
addition, the data are ordered first by bin number and then by time within each bin. The surface 
elevations on this data set have the orbit adjustment applied where it was available. If the orbit 
adjustment was not available, (indicated by the orbit adjustment value for that record being set 
to -9999) then the surface elevation contains the value calculated from the unadjusted orbit. The 
slope correction has not been applied to any of the surface elevations. 

The data base is designed to be used on a direct-access device, so that data from one or 
several bins may be accessed without the need to read all the records prior to the location 
desired. This is achieved by dividing the data base into three sections. 

The first section of the data base, a header, consists of one logical record and gives a 
summary of its configuration: the locations of the corners of the data base, the number of 
latitudinal rows, the width in degrees of each of these rows, and the number of longitudinal 
divisions in each row. These pieces of information give the layout of the data base, as depicted 
in Figure 10. Information pertaining to the size of the data base, the starting record of the bin 
directory, and the corrections applied to the data are also contained in this header. 
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Figure 10. Seasat Greenland Geo-referenced Data Base Configuration 
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Following the header are the altimetry data ordered by bin number and, within each bin, 
by time. The altimetry data are subdivided into two groups for each bin which contains data. 
The first subgroup consists of one logical record which indicates the number of data points 
contained in the bin. The second subgroup consists of the actual altimetry data (position, rev 
number, surface height, orbit adjustment and slope correction), with each record corresponding 
to a data point. 

The final section is a bin directory which follows the altimetry data. The bin directory 
starts at the logical record indicated in the data base header. The directory contains an entry for 
each bin. and starting with the first bin, indicates the record number in the data base (not 
including the header record) at which the start of the data for a particular bin may be found. 

Bins which contain no data have a zero entered in the directory. Table 7 summarizes the 
structure of the data base in greater detail. 

One use of the data base is to assist the gridding program (Section 5.0) in locating and 
accessing all data contained within a specified radius of a grid location. In addition, the data 
base may be used to locate data within any desired area. The following example demonstrates 
how this may be done. The limits of a desired area are used in conjunction with the header 
information to determine exactly which bin numbers contain the data. Using the southernmost 
latitude of the desired area, along with the width of the latitude rows, establishes the 
southernmost row which contains the data. Longitudinal limits of the desired area are then 
checked In conjunction with the size and location of the longitudinal divisions in that row. When 
the longitude limit of the desired area for that latitudinal group is exceeded, the process starts 
again with the next latitude row to the north. These steps are repeated until the northernmost 
boundary limit of the desired area is reached. 

Equipped with the bin numbers which contain the data, the directory, which gives the 
logical record on the direct-access disk at which each bin begins, is read. If the directory value 
for the bin is non-zero, this logical record is then read to determine the number of records which 
follow and are contained in the same bin. The subsequent data is then read for each bin. 
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SECTION 5.0 
GRIDS 


The uneven distribution of Seasat data presents problems when attempting to create 
computer generated contours. An intermediate step is useful which fits the data to nodes of a 
regular grid. Data local to each grid point are fit with a biquadratic or bilinear surface to 
determine the surface height at the grid point. This procedure is referred to as gridding the data. 
Grids are generated using the corrected and adjusted data in the geographical data base. 

5. 1 POLAR STEREOGRAPHIC PROJECTION 

Grids of the Greenland data are generated in a tangent polar stereographic projection 
where the plane of projection is located at the geographic North Pole (the projection latitude) 
and is normal to the earth's axis. This projection is conformal which results in equality of scale 
about a point. Figure 11a depicts the concept behind this type or projection. A straight line is 
drawn from the South Pole (pole of projection), through a point on the earth’s surface, 9, to the 
projection plane which is tangent to the North Pole. The projection plane is in turn divided into 
square grids from the pole to the Equator with the North Pole at the center. Three projection 
parameters define the size and the orientation of the plane and the grid size: 

S - a conversion factor from half-inch grids at the projection latitude to the desired 
grid size; 

4> p - the minimum latitude extent of the map perimeter for the projection latitude located 
at the North Pole; maximum latitude extent for the projection latitude located at the 
South Pole; 

G - the Greenwich orientation in degrees 

In the case of Greenland, where 20 kilometer grid cells were decided as being optimum for the 
data distribution, values of S= 1.65, <f> p = 50°, and G=45° were chosen. 

These three parameters are sufficient to define a grid of the northern hemisphere, from 
the North Pole to 50° north latitude where the number of cells of desired size from the pole to 
the equator may be represented by: 
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Figures 1 la and 1 lb. Polar Stereographic Projection of Point Q with Latitude <f> 
and Longitude A onto Plane with Map Perimeter 5CP 





D 


2R 


(5.1) 


S x 10 6 


where R is the radius of the earth measured in one half-inch grid cells and was chosen to be 
consistent with polar stereographic projections described in other documentation. 


The integer number of grids of desired size from the pole to the map perimeter is. 


90 - |<f> I 
N = D x tan IL21 

2 


(5.2) 


The grid, defined by I and J axes, with the origin in the upper left corner (see Figure 
lib), represents the coordinate of the North Pole as: 


Jp = N + 1 . 


Any point with latitude 4> and longitude A, which is located in the northern hemisphere 
north of is positioned at the following I, J coordinate: 


I = INT [d x A x cos(X) + Ip + 0.5] 

j = INT [d x sin(X) + Jp + 0.5) 


(5.4) 


where 


d 


is D x tan 


90 - fej 
2 


X is A + G 


A is 4- 1 if 4> p ^ 0 
A is -1 if <fi p < 0 . 
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5.2 GRIDDING PROCEDURE 


The surface height at each grid point location is calculated by fitting the surrounding 
data to the following biquadratic surface modeling function: 


h,j (A ,4>) 


Ci + C 2 ( A - A f ) + 
*J 'J capmin 


^3 [<f> ~ <f>jl 

lJ (capmin)(cos(£j) 


where 


h 


ij 


+ c 4 (A_-.A,) - 4>A + c 5 (A - A,) 2 

1J (capmin) (capmin)(cos^>j) capmin 2 

+ C 6 {<t> - 

•* (cos (capmin 7 ) 


(5.5) 


= value of the surface elevation function for the ij grid point as evaluated at 
the location (A, <fi); 


° 'ij _ ° 6 ij = numerical| y determined coefficients or the biquadratic function for grid point 
ij; and 

= longitude and latitude of the ij grid point in deg. 
capmin = minimum cap size in deg longitude. 


A weighted least-squares method is used to solve for the coefficients Cj - C 6 at each grid 
point ij. The weighting is invoked to prevent the obliteration of the local surface details by the 
smoothing process, and to lend greater importance to the data closest to the grid point location. 
The form of the weighting function is 



Ou D,. 


(5.6) 


where 



weight of the k th data point used in determining the coefficients of the surface 
function for the ij grid location; 

observation standard deviation of the k th data point; 


N = power of inverse distance weighting; and 

D k ^ = the distance from the k th data point to location ij, 
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where 


D kij = {l(A k _ A i } cos ^k 1 + (< ^k " } 


2i 1/2 


The observation standard deviation was assigned a value of 1.0 m. The power of inverse 
distance weighting was assigned a value of 2.0 m. The formula used for the least-squares 
minimization in matrix notation is 


p !j w ij p u c ij = p ij w u H u 


(5.7) 


or 


C u = [ p ij W u p .j] _1 p !j W,J H,j 


(5.8) 


where 


H U = 


IS 


the observational data set used in determination of grid point ij; 


p u = 


ah, 


ac, 


2h r 

ac. 


ah, 


ac. 


ah, 


ac. 


2h 


HL 


ac. 


is the matrix of observational partial derivatives; 


is the set of coefficients for grid point; 
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o 


w = 

w Ij 


w, 




u 


is the observation weighting matrix. 


A solution exists for Equation (5.8) if the determinant of the normal matrix WyC l 

is positive. However, poor data distribution can cause ill-conditioned matrices yielding solutions 
that vary considerably from the expected results. One needs to be able to recognize when 
numerical problems occur to assure reasonable solutions. To this end the singular value 
decomposition (SVD) method is used to solve the matrix equation. The results of the SVD process 
give an indication of the stability of the equations and therefore whether a unique stable solution 
exists (Forsythe, Malcolm, and Moler, 1977). When the normal matrix By is used as input to SVD, 
three output matrices are calculated: 2 , U, and V. 2 is a diagonal matrix, such that 


2 = 


where the o'a are referred to as the singular values of B. The matrices U and V are used to 
transform the equations 

Bc=y 

into an equivalent diagonal set of equations 
2 c = y . 


In principle, if none of the a’s are zero the transformed equations could be solved using 


C, 


y ( 


a 


l 


In practice, when any of the a's are small, numerical instability can result, giving unreasonable 
answers. The key to using SVD is to set a tolerance r which reflects the accuracy of the data 
and the arithmetic used. If any a’s are less than r times the largest a then those corresponding 
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1T-S are not uniquely defined and unreasonable results can occur. When problems occur, steps 
must be taken to provide more information to evaluate the surface function. 


Once x is chosen, then 2, U, and V are used in the following manner to calculate each 
coefficient 


m 

S = 2 U(|,i) Yj 

j=i 


for all j where > t 


c . 


J A V(i,k) . 
k= 1 a k 


In this study the value of r used was .001 m. SVD is then used to determine when there 
are sufficient data to provide a unique solution to the surface modeling function. When a unique 
solution cannot be found more data are added and the function is reevaluated. At each grid 
location ij, data within the circular area defined by radius R from the grid location are used in 
the solution. Four different values for R are used: 33 km, 55 km, 88 km, and 132 km. Initially 
the smallest value of R is used and if a solution cannot be found then R is increased. If the 
biquadratic solution at the maximum value of R is unsatisfactory according to the SVD criterion, 
then the function (Equation 5.5) is reduced to a bilinear function by setting coefficients C4 
through C6 to zero. If a valid solution still cannot be found, then the grid value is considered 
undefined and set to -1000.0. 

Individual data point removal is also invoked during the gridding process. After finding 
a valid solution at location ij, the weighted rms of the residuals of the data with respect to the 
surface is calculated using 

2 

M Res M) 

2 _2 N 

k=l °k D k|j 

2 

M l 

2 a 2 N 

k=l °k% 


RMS 


WT = 
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where 


Res kij ~ h k " h kij 


h kij ~ height at location of measurement k evaluated using the surface function 

for grid location ij. 

The following inequality is then evaluated for each data point used in the solution. 

< E mult * RMS WT (5.10) 

a N J 

°k D 

A value of the editing multiplier (E mult ) equal to 3.5 is used and all data points that do not 
satisfy the inequality are deleted. When any data points are deleted the surface function is 
reevaluated using the remaining data. A minimum of 10 data points are required to solve for 
the function. 

The standard deviation associated with the grid height, a Gt j t is then calculated to determine 
how well the grid represents the data. 

a GIJ = RMS WTjj * (VllU ) 1 ^ 2 ’ 


where 


v u = 


B -1 P W 
tj w u 


°i 


“m 


L 0 


[B7j PTj Wy] 


Grid points that have a large value of a Q do not represent the data as well as those that 


have smaller a G 's. 


The format of the grid record is described in Table 8. The location. coefTlcients 
number of points used and other pertinent parameters are output for each grid point location. 
The user can utilize these parameters to decide the accuracy of the individual grid values. 
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Table 1. Ice Data Record Description 


General Characteristics: 

Record Format - variable 

Record Size (bytes) - 164 + 4 for IBM record control word 
Blocksize (bytes) - 31920 + 4 for IBM block control word 

The first seven records of the IDR data set are 80 bytes long and contain a brief description of 
the contents of the file. The remaining records follow the 164-byte format. 

HEADER RECORDS 


FORTRAN 

Variable 

Bvtes Type Description 

1_30 Ai Brief description of file contents. (Comprises first seven records 

only) 


DATA RECORDS 


FORTRAN 

Variable 

Bvtes Type Description 

!_4 1*4 Satellite ID - This is the international satellite designation nnpppqq 

where: 

nn - last two digits of the year of launch (e.g., 1974 74, 1969 
69). 

ppp - order of launch. 

Example: 

The 25th vehicle launch in a given year is designated with 
ppp = 025. 

qq - component identifier (e.g., 
component a -* 01, 
component € -► 12, etc.). 

5-6 1*2 Measurement type 

40-44 Altimeter height 

40 = Long pulse (GEOS data) 

41 = Short pulse (GEOS data) 

43 = Seasat altimetry 
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Table 1. Ice Data Record Description (Cont.) 


Bytes 

7-8 


9-12 

13-16 


FORTRAN 

Variable 

-.Type 

1*2 


1*4 

1*4 


Description 

Time system indicator (nm) 

n-value Description Tracking data times 

0 Ground received time 

1 Satellite transponder/reflector time 

2 Ground transmitted time 

Altimeter data times 

1 Transmitter time 

2 Ground bounce time 

3 Receiver time 

m -value Description 

0 UT-0 

1 UT-1 

2 UT-2 

3 UTC 

4 A. 1 

5 A.3 (A.T. B.I.H.) 

6 A-S (Smithsonian) 

Station Number (0 indicates altimeter) 

Altimeter measurement status word 

The status word consists of bit switches packed into a single 32- 
bit word. The rightmost bit (bit 31) is of lowest order and the 
leftmost bit (bit 0) is of highest order. 


0 


31 


The status bits are configured as follows: 

Bits Value Description 
0 Unused 


1 

2 

3 


Format indicator for measurement types 40-4 1 
20 obs/frame 
32 obs/frame 
320 obs/frame 



Table 1. Ice Data Record Description (Cont.) 


FORTRAN 

Variable 

Bytes Type Description 

(13-16 Cont.) Bits Value 

3 

0 

1 

4 

5-6 

0 

3 

7 

8 

0 

1 

9 

0 

1 

10-11 

0 

1 

12 

0 

1 

13 

0 

1 


Description 

Net instrument corrections indicator 
Instrument corrections applied to observation 
Instrument corrections not applied 

Unused 

Speed of light indicator 
2.997925xl0 8 meters/sec 
2.99792458x1 0 8 meters/sec 

Unused 

Solid tide indicator 

Solid tide not on data record 

Solid tide on data record 

Ocean tide indicator 

Ocean tides not included in total tides 

Ocean tides included in total tides 

Tropospheric correction indicator 

Total tropospheric correction not on data 

record 

Total tropospheric correction on data record 

Ionospheric correction indicator 
Ionospheric correction not on data record 
Ionospheric correction on data record 

Atmospheric corrections indicator 
Ionospheric and tropospheric corrections 
applied to observation if found on data record 
Ionospheric and tropospheric corrections not 
applied to observation if found on data record 
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Table 1. Ice Data Record Description (Cont.) 


FORTRAN 

Variable 

Bytes Type Description 

(13-16 Cont.) Bits Value 


14 

0 

1 

15 

0 

1 


16-20 

21 

0 

1 

22-27 

28 

0 

1 

29 

0 

1 


30 

0 

1 


31 

0 

1 


Description 
Total tide indicator 

Solid and ocean tides removed from 
observation if found on data record 
Observation includes solid and ocean tides 

Center of gravity indicator 

Center of gravity correction applied to 

observation 

Center of gravity correction not applied to 
observation 

Unused 

Altimeter mode (GEOS only) 

Global track mode 
Intensive track mode 

Unused 

Location indicator 
Over water 
Over land 

Orbit adjustment indicator 

Orbit adjustment has been applied to 

observation 

Orbit adjustment has not been applied to 
observation 

Slope correction indicator 

Slope correction has been applied to 

observation 

Slope correction has not been applied to 
observation 

Retracking correction indicator 

Retracking correction has been applied to 

observation 

Retracking correction has not been applied 
to observation 
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Table 1. Ice Data Record Description (Cont.) 


Bvtes 

17-20 

21-28 

29-36 

37-40 

41-44 

45-48 

49-52 

53-56 

57-60 

61-64 

65-68 

69-72 

73-76 

77-78 


FORTRAN 

Variable 

Type 

1*4 

R*8 

R*8 

R*4 

R*4 

R*4 

R*4 

R*4 

R*4 

R*4 

R*4 

1*4 

1*4 

1*2 


Description 

Modified Julian Date (MJD) of observation 
Julian Date = MJD + 2400000.5 

Fraction of day past midnight (GMT) 

Altimeter range measurement in meters 

Satellite latitude in degrees 

Satellite east longitude in degrees 

Measurement standard deviation in meters 

Center of gravity correction in meters 

Tropospheric refraction correction in meters 

Ionospheric refraction correction in meters 

GEM10-B geoid height above reference ellipsoid in meters 

Total tide height above reference ellipsoid in cm. 

Rev number 

Surface height with respect to ellipsoid in cm. 

Surface height status word 

-I 

15 

Description 
Unused 

Slope correction applied 
Slope correction not applied 

Orbit adjustment applied 
Orbit adjustment not applied 

Solid tides removed 
Solid tides not removed 

Retracking correction applied 
Retracking correction not applied 

Center of gravity bias applied 
Center of gravity bias not applied 


0 


Bits Value 

0-8 0 

9 1 

0 

10 1 

0 


11 1 
0 


12 1 
0 


13 1 

0 
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Table 1. Ice Data Record Description (Cont.) 


Bytes 

(77-78 


79-80 

81-82 

83-84 

85-86 

87-88 

89-90 

91-92 


93-96 

97-100 


FORTRAN 

Variable 

Type 

Cont.) 


Bits 

Value 

Description 

14 

1 

Tropospheric correction applied 


0 

Tropospheric correction not applied 

15 

1 

Ionospheric correction applied 


0 

Ionospheric correction not applied 


**2 Significant wave height (H 1/3) in cm. 

J * 2 Automatic Gain Control (AGC) in dB 

1*2 Solid tides in cm. 

l * 2 Tangent of along-track slope (x 10 5 ) 

I* 2 Tangent of cross-track slope (x 10 5 ) 

J * 2 Size of window used in obtaining along-track slope in meters 

Along-track and cross-track slope correction word. If all bits are 
zero, then slopes for slope correction were not able to be 
computed. 


1 

0 


— i 

15 

Bits 

Value 

Description 

0-9 


Unused 

10 

1 

Along-track slope set to the maximum value 
of .8 degree during iterative procedure. 

11 

1 

Cross-track slope set to the maximum value 
of .8 degree. 

12 

1 

Along-track slope set to .8 degree after final 
iteration. 

13 

1 

Window was extended to 20 km with no point 
found; reference grid used to calculate along- 
track slope. 

14 

1 

Window had to be extended in both directions 
to determine along-track slope, but it is less 
than 20 km. 

15 

1 

Two adjacent points were found and used to 
determine along-track slope. 


Orbit adjustment to 84,306 ocean surface in meters 
RMS of orbit adjustment fit in meters 
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Table 1. Ice Data Record Description (Cont.) 


Bvtes 

FORTRAN 

Variable 

Type 


101-104 

R*4 

RMS of filtered fit in counts 

105-108 

R*4 

Timing bias in seconds 

109-144 

R*4 

Retracking parameters 
>3(1) ->3(9) 

145-148 

R*4 

Attitude information from SDR 

149-152 

RM 

Correction to surface height if using leading edge of leading edge 
in meters 

153-156 

1*4 

Geographical data base bin number 

157-158 

1*2 

Standard deviation of 1st leading edge position in gates 

159-160 

1*2 

Standard deviation of 2nd leading edge position in gates 

161-162 

1*2 

Retracking status word 


0 


15 


Bits Value Description 
0 Unused 


2 

3 

4 

5 


0 Gains and offsets were not applied to waveform 

counts in plots and in determining ft 

parameters 

1 Gains and offsets were applied to waveform 

counts in plots and in determining /3 

parameters 

0 Specular test not performed or waveform not 
specularly shaped 

1 Waveform determined to be specularly shaped 

0 Status flag from SDR less than or equal to 
one 

1 Status flag from SDR greater than one 

0 Waveform not specularly retracked 

1 Waveform specularly retracked 

0 Gains and offsets not applied to waveform 
count values on WDR's 

1 Gains and offsets applied to waveform count 
values on WDR’s 
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Table 1. Ice Data Record Description (Cont.) 


FORTRAN 

Variable 

Bytes Type Description 


(161-162 Cont.) 


applies to 
water data 


163-164 1*2 


Bits 

Value 

Description 

6 

0 

For double waveforms the retracking correction 
is not calculated from a weighted average of 
the two leading edges. 


1 

For double waveforms the retracking correction 
is calculated from a weighted average of the 
two leading edges. 

7 

0 

No problem with leading edge definition of 
waveform 


1 

Waveform not defined well enough to filter, 
no leading edges or too many leading edges 

8 

0 

No problem retracking 


1 

Problem retracking 

9 

0 

Timing bias was not applied to time tag 


1 

Timing bias applied to time tag 

10 

0 

Waveform not retracked 


1 

Waveform retracked 

11 

0 

Whole edge retracked 


1 

Leading edge retracked 

12 

0 

Ht correction not applied due to h 


1 

Ht correction applied due to h 

13 

0 

Attitude seastate correction not applied to h 


1 

Attitude seastate correction applied to h 

14-15 

0 

Tracking mode 1 


1 

Tracking mode 2 


2 

Tracking mode 3 


3 

Tracking mode 4 


Version number of retracking program that converted the data 
from SDR to IDR format 

n l n 2 n 3 n 4 n 

5 

P 

^P 

to 

11 

year of version 

n 3 n 4 = 

month of version 

n 5 

point no. of version 
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3522 3213 3112 

574 71.87 D 72.064 324.774 71.898 315.464 33 4262 4258 4140 4139 4039 


Table 2. Seasat IDR Greenland Catalog (Cont.) 
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Table 4. Corrections To Seasat Ice Data Records 


CORRECTION 

ADJUSTMENT 

VALUE OR 
RANGE 

MANNER IN WHICH APPLIED 

SECTION 
IN WHICH 
DOCUMENTED 

TIME 

SURFACE 

HEIGHT 

RETRACKING CORRECTION 
accounts for lag in 
tracker response 

— 1 5m < Ah 15m 

N/A 

(-) 

2.1 

TIME BIAS accounts for 
track mode correction 

-7.9451 x 10~ 2 s 

(+) 

N/A 

2.2.1 

SIGNAL TRAVEL TIME 
CORRECTION 

-2.67 x 1CT 3 s 

(+) 

N/A 

2.2.1 

CENTER OF GRAVITY 
OFFSET adjusts 
measurement to s/c 
center of mass 

~ 6.04 m 

N/A 

(-) 

2.2.2 

IONOSPHERIC REFRACTION 
CORRECTION accounts 
for signal delay 

~ 2-3 cm 

N/A 

i+) 

2.3.1 

TROPOSPHERIC REFRAC- 
TION CORRECTION 
accounts for signal 
delay 

~ 1. 5-2.5 m 

N/A 

(+) 

2.3.2 

SOLID TIDE removal 

~ 2-10 cm. 

N/A 

h 

2.4 

ORBIT ADJUSTMENT 
reduces orbit error 
and references the 
data to a mean ocean 
surface 

3m< AH orb < 3m 

N/A 

(-) 

2.5 

SLOPE CORRECTION 
accounts for signal 
being returned from 
closest point within 
satellite footprint 

Om<AH SLOPE <80m 

N/A 

(-) 

2.6 

i 
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Table 5. Waveform Data Record Description 


General Characteristics: 

Record Format 
Record Size (bytes) 
Blocksize (bytes) 


FORTRAN 

Variable 

Rvtes Tvpe 

1-8 

R*8 

9-16 

R*8 

17-20 

R*4 

21-24 

R*4 

25-28 

R*4 

29-32 

R*4 

33-36 

1*4 

37-38 

1*2 

39-40 

1*2 

41-166 

1*2 

167-168 

1*2 


variable 

170 + 4 for IBM record control word 
31842 + 4 for IBM block control word 


Description 


Fraction of day past midnight from sensor data record 
Altimeter range measurement in meters from sensor data record 
Satellite latitude in degrees from sensor data record 
Satellite east longitude in degrees from sensor data record 
Altitude error Ah in meters 


Altitude rate error Ah in meters/sec 

Modified Julian Date of observation from sensor data record 
Significant wave height (H 1/3) in cm. 

Automatic Gain Control (AGC) in dB 
Waveform counts 

Word indicating original data flags 


Bits 

Value 

Description 

0-10 


Unused 

11 

1 

Not in track mode 

12 

1 

Chirp/cw 

13 

1 

Altimeter error status 

14 

1 

Reacquisition 

15 

1 

Acq/Trk 
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Table 5. Waveform Data Record Description (Cont.) 


FORTRAN 

Variable 

Bytes Type Description 

169-170 1*2 Retracking status word 


0 


15 


Bits Value Description 
0 Unused 


1 


2 


3 

4 

5 


6 


0 Gains and offsets were not applied to waveform 

counts in plots and in determining (3 

parameters 

1 Gains and offsets were applied to waveform 

counts in plots and in determining /3 

parameters 

0 Specular test not performed or waveform not 
specularly shaped 

1 Waveform determined to be specularly shaped 

0 Status flag from SDR less than or equal to one 

1 Status flag from SDR greater than one 

0 Waveform not specularly retracked 

1 Waveform specularly retracked 

0 Gains and offsets not applied to waveform 
count values on WDR’s 

1 Gains and offsets applied to waveform count 
values on WDR’s 

0 For double waveforms the retracking correction 
is not calculated from a weighted average of 
the two leading edges 

1 For double waveforms the re tracking correction 
is calculated from a weighted average of the 
two leading edges 
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Table 5. Waveform Data Record Description 


FORTRAN 

Variable 

Bytes Type Description 

169-170 Bits Value 

(cont.) 

7 0 
1 


8 0 

1 

9 0 

1 

10 0 

1 

11 0 

1 

12 0 

applies to 1 

water data 

13 0 

1 

14-15 0 

1 
2 
3 


Description 

No problem with leading edge definition of m 
Waveform not defined well enough to filter, no 
leading edges or too many leading edges 

No problem retracking 
Problem retracking 

Timing bias was not applied to time tag 
Timing bias applied to time tag 

Waveform not retracked 
Waveform retracked 

Whole edge retracked 
Leading edge retracked 

Ht correction not applied due to h 
Ht correction applied due to h 

Attitude seastate correction not applied to h 
Attitude seastate correction applied to h 

Tracking mode 1 
Tracking mode 2 
Tracking mode 3 
Tracking mode 4 
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Table 6. Seasat Greenland Geographical Data Base 


A> 



rs r-\ 

rv 

A> AN A' A(, 

r-\ r\ 

rv rv 


r, 


NO 

i — 4 

o 

LT\ CNJ 

NO 

KD CM NX CD 

O' NO 

ON NO 

rH 

O' hi 

CM 

H 


pH 

rH r— 4 


IAK) rHCM 

CM 


CM 

rH 

w 


w 

vy w 

<Ny 

N>UN> W 

n> 

w w 


Ny w 


IA 

vX 

oo 

00 rH 

CM 

00 vT iA 00 

0O rH 

LA O' 

CM 

LA 00 

iA 

A* 

On 

OO 

00 IA 

<NI 

A- ON NX r-H 

oo VO 

NO A- 

CM 

NX to 

A- 

(NJ 

*0 

00 

00 NO 

CM 

oo fo OO fo 

0O NO 

0O rH 

CM 

CO oo 

CM 

r-H 

H 


rH 

rH 

rH rH 

rH 

rH 

rH 

rH 


r> 

AN 

AN 

As As 

an 

rsAA/NA 

AN As 

r> r> 


pH r-v 

an 

00 

O 

CM 

CM vo 

pH 

IA rH NO IA OO 

CM A- 

00 CM 

tA 

A CM 

NO 

CM 

rH 

CM 

CM 


NX 

NO 

IA 

CM 

CM 

Ny 

sy 

w 

w 

W 

vy vy Ny sy w 

*y sy 

w w 

w 

w w 

sy 

CM 

rH 

iA 

00 IA 

O' 

IAIA rH A- iA 

o 00 

O' CO 

O' 

O LA 

(M 

hi 

IA 

vX 

A- A- 

A* 

fOXIA AN 

CM rH 

IA A 

A 

CM vX 

NO 

CM 

NO 

00 

00 CM 

rH 

oo OO NO NO CM 

A- NO 

A 00 

pH 

A 00 

CM 

rH 

rH 


rH 

rH 

r-H ri 

rH 


rH 

*H 




r\ 

r. r\ 

an 


AS 

pH rH 

H 

HHHHH 

H H, 


AAArHryNA/iN 


/a, 


O 

A O 

o 

pH 

IA 

rH rH 

A 

CM On IA NO pH 

CO CM 


OOpHnXpHONOCMCM 


tO 

J— 



vX 

pH 


rH 



NO CM IA nX 

CM CM 


CM rH NO rH 


a. 


n> 

sy s_y 

sy 

w 

w 

w sy 

sy 

Ny w ^y ^y sy 

Ny Ny 


Ny <u> Ny w <^y Ny Ny Ny sy 




NX 

OO Ov 

oo 

OO 

NX 

A iA 

A 

CM vX OO vX IA 

A IA 


vXCMihCMAvX'XOOv 


QC 


o 

oo oo 

CD 

CD 

NX 

A vo 

0O 

O vX O nX so 

00 A 


NOOIAONOvXvXOCO 


UJ 


NO 

00 rH 

NO 

NO 

NO 

VO CM 

NO 

00 vo NO nO CM 

NO CM 


NOOONDOO'X'ONO'XrH 


CO 

s: 

3 


rH 

pH 

rH 

pH 


pH 


rH pH 

rH 


»H pH rH pH 


z 

w 

an 

AS 

As As 

r. 

yN 

r> pH 

as r> 

hhh 

hh y>y> hh 

rn r> 


A A A A AAAAAA 

AN 


NO 

NO 

NO A 

vX 

NO 

CM rH 

rH pH 

CM CM CM 

NO NO no CM A vO 

NX O' 

nX 

IT| vo ro K) so (\J VO H O fsl 

NO 

> 

UJ 



pH 


pH 

rH CM 

NX 


rH rH CM NO 



pH pH pH iH 


a: 

sy 

sy 

sy Ny 

>sy 


sy sy 

sy sy 

w sy w 

Ny sy Ny Ny W ^y 

<wy w 

w 

*ysyNyNy«Nywsy^yNyNy 



NO 

r—i 

CM 00 

OO 

uo 

IA O 

NO CM 

A vXO 

IA»X.HlAOCM 

vo CM 

NO 

OOAOOAvXrHOAAA 

NO 


A 

iA 

C7V 0O 

A 

NO 

VX O 

vX NO 

A O' CM 

NO NO O NO O NO 

A NO 

vO 

IA A O A ON O O nX pH nX 

pH 


A 

NO 

A OO 

OO 

CM 

0 O vX 

vX CM 

NO NO A 

00 NO NO CM vX CM 

NT CM 

A 

IAnONOnONOvOvXnXOnX 

A 



pH 



rH 


rH 

pH 

*H pH 

pH 


H pH r-i i—i r—i 



as 

an rn 

pA rn 

AN 

/A 

AA A A A AA AA 

a an a a a a a rn rn a an an a a a a a a a a a a a 


vX 

IA O' 

IA NO 

LA 

CM 

CMOVpHnXOVNOOCMpH 
NO rH 

'OACMOvUOvOvOrHNOvOUOpHrHrHK>rHONOOmOVAOO 

pH pH CM r—i r—i pH rH 

Ch 

s_y 

W Ny 

Ny N_y 

Ny 

Ny 

Ny w Ny Ny w Ny Ny Ny Ny 

Ny Ny Ny Ny Ny N_y Ny Ny Ny Ny Ny Ny sy Ny Ny Ny w Ny Ny w n> vy 


O' 

ON CO 

rH A 

O' 

O' 

vX vO CM O' Ov A NO H no 

CvJpHOCMO'AOfOO'OhOO'O'pHtOLJOrHOO'OLOCMApH 

CM 

CM 

CM CD 

OV 0O 

00 

CM 

'XiANOcOOONO'XiAA 

OOOCMcOrHO'vXCOO'vOrHUOAiA'OOCDAAK)oOvo 

rn 

iA 

LA NO 
pH 

LA vO 

pH 

IA 

VOr-tCMrHrH'TvXNOvX 
rH rH rH rH 

eOvOvXCMpHONXvXpHNXAUOACMLTlCMvo^vXCMCMvOfO 

pH pH pH pH pH r-i rH rH pH 

vX 

rn rH 

A A A A A A A 

AAA A 

rn rn 

rn a rn a rn as fa a rn a rn a 

rn an an rn rn a an a rH an an a an an a an an an a an an an an an 

AN 


(M<MrHNOrHCMCMNoosr^o<MN , >'xcMvocMrH'XNOvor'*oorHNOr-icoCMiricoinvxco^'(Mov\oNOrHCMiri'OrH'XO'XoviriNocMov.HCM 

r-» NO r-i r-i CM rH i— I pH CM vX *— I •— I r- I i— I r- I <M CM CM (\J 


<* Ny v> vy n 


o>Lnc\i*— »ir>sooc\ioK>«— <\3'v©r>-cr'CT' , eo»— <•— ivoowofNwNxovirtcMp— ioovovoovNocoi s «.CMr ,,, '.F*'*iciinr- tvovoCMovtocMCMovinoosoO' 
CM»r^voNOvxcMo' 0 'XsoNONXNOoo(Muooo'vTooin»“« 0 'flOvoNOoMri'XvXr^irioo»ONOrH'Xr'.f'-ov'X'X<McO'ONOK)oor^-*Hfv.(\j 

iOCJLnNOCOrHlOr'.-NXvTNOvOp— <vTpHtnirivOinr-ip--*r-tON , '>iHCMCMiJOiniHpHpHrHOOvXCMOvXCMCM4/'lpHrH<MrH(MCMCMp— I CM NO NX vT 


OC 

UJ 


QC 

o 

o 


oooooooooo oo 

oooooooooo oo 


o o o o o o 


o 

o 


o o o o o o o 

O CD CD O O O O 


O O O O O 
CD O O CD CD 


O O CD O O 

O O O O CD 


3 

'XLn'OCMvXtnvOrHCMHO 

vX uo 

NO r^- o rH CM NO 

vX 

LO 

NO o rH CM NO 

vT 

to 

pH pH pH r—i r—i pH t — 1 r—i r—i rH 

ri pH 

pH r—i pH r—i r—i i — i 

pH 

H 

rH »H pH r—i rH 

pH 

CD 

hONONONO^OfONOKONOrO 

NO NO 

NO NO NO NO NO NO 

NO 

NO 

NO NO NO NO Nl 

NO 

z 

a 

oooooooooo 

O CD 

oooooo 

o 

o 

O O O O O 

O 

— » 
i 

O'O'O'vX'XvX'TO' O' O' 

O' O' 

O' O' vX vX vX O" 

vX 

vX 

vX Ov O' O' O' 

O' 

»- 

O'O'O'OOOOOOO 

O O 

O O pH rH rH pH 

rH 

pH 

pH rH pH rH pH 

r—i 

< 

Lmnmsovo'O'ONO no vo 

vO vO 

NO vo 'O 'O vO vO 

VO 

VO 

vo NO NO NO vO 

vO 


m 

\OSOHCJ 

NO 

vX iO 

SO 

CO O' 

pH 

rH pH rH rH rH 



pH 

CD O 

NO 

NO NO NO tO NO 

NO 

NO NO 

NO 

NO NO 

CD 

O O CD O CD 

CD 

o o 

CD 

a o 

O' 

O' O' VT vX VX 

vX 

'X'X 

vX 

oo 00 

pH 

pH rH CM CM CM 

CM 

CM CnI 

CM 

CM CM 

SO 

vO vo vQ 'O vo 

'O 

vO SO 

vO 

'O vo 


00 


CL 


UJ 


flq 

£ 


Ma'HNOHNvoro 
rH CM 


> 3 “ NT 
mfM 


cm 

NO 


fh vo CM p-» oo r". vx 

CM CM O fv. 


rH >H CM vT IA sX 00 

NT r-i (NJ NO H h- 

NO CNJ 


lO 

OO 


nONnCON 
H IA >H 


h- o 
rn*o 


IA 

r-'- 


or 

UJ 

« 

3 


IA NO »0 IA VO A- CNJ M NT IA so 

rH rH i-H IA LA 1A IA On On On O' On 


A* OO H (NJ *0 nX IA 
ON O' N*> eo to hi K) 
H ft pH H H 


NO 

NO 


fs. rH CM NO 'X 
HNM--N 


m 

so 

fs. OO rH CM NO 

vX 

l n vo 

1^. 

O' CD 

r- 

Cs- 

r*** S pH pH pH 

pH 

pH pH 

r—i 

■r m 

p-t 

iH 

iH rH CM CM CM 

CM 

CM CM 

CM 

CM CM 


92 


r— I | — K>i-t CO VO 

X rHCM CM 

w w ■w w *>-> ^ 

Ov rH <M CM in CM 

in in o o x m 

r>. m oo co co cm 


AN AN 

AN AN 

AN AN AN AN 

AN 

r*\ 

O' CM 

rH rH 

m m vo m 

m 

O 


rH CM 

m rH 

rH 

H 

NA W 

N/ VA 

W NA NA NA 

NA 


O' rH 

r- oo 

o o m O' 

rH 

in 

m m 

o 

CM CM vo eo 

vO 


fN. m 

so m 

fN. jn. 1 "- rH 

m 


rH 

rH 

rH 

rH 

H 


AN AN 

fN. NO 

x m 

55) 

(IT 

5) 

AN AN 

mcM 

AN AN AN AN 

ah on eo Ov 

rH m m cm 

An 

vo 

11) 

16) 

7 59 C 
720 ( 

CM 

o 

oo 

1318 ( 

NA 

r- 

X 

X 

rH 

1 1 7 5 C 
1351 C 

o o mm 
CM CM vO so 

va 

rH 

VO 

m 

rH 

na 

X 

O' 

m 

rH 

1 175( 

AN AN 

X VO 
X 

(61 

iSL 

15) 

5 ) 

AN AN AN 

m m in 

AV AN AN A AN 

m rH m fn. x 

NO CM 

(I 

10) 

AN AN 

X CO 
rH 

o «0 
cm m 
r- m 

NA N/ 

mcM 
vo o 
r- co 

1 275 ( 

rH 

VO 

m 

rH 

NA NA NA 

X 00 co 
crv r- o 

m co m 

rH rH 

W WA NA SA W 

VO vo ON OV vO 

rHrHinirtO 

r^- r- r- co 

1 318 C 

1351 ( 

w NV 

rH X 

vo r-» 
CM 00 

rH 


AN AN ^v I^V ^ ^ ^ 

(M M) tON'OONvJ 

r- 1 K)'OHH 


00 m OOHVJNOCO 

LTl SO lAOvTHO'H 

in cm in \o vo o x no 


rH NO CM 

O VO m 

vo cm 


an an an an an 

vo vo oo oo r- 
N»o*o 


VO sOOv CM CM 
r*^ r^H voo 
vx x m moo 


oo oo 

O CM CM 

m r- h- 


rv an an rv (TV ^ an ^ 
mo vo cm sr r- o 

CM HH r~» r-4 


AN AN AN 

(v.Nj<r vr cm h- m 

rH rH rH CM 


O' m m x x 

m m rH rH 


m n*. n") cm vx »-h vo sr 

hmio^h 


m cm vo O' oo r ^- cm 

rH CM f-- rH rH OO X m 

m cm x m m oo x cm 


K) VO CM CM m VO vo 
r- r'- r- n n 

nO X X X CM X X 


ton 0 o m 
mo »^n 

m VO rH X CM 


X m rH CM CMm© 0 — « 

r^. r-. vo |s.NHifto 

oo vo cm vx vx m m vo 


A' /-V AN AN AAA 
SNf CO'ONN'OmX 
rH X m K> rH 


rH-OV CM vo VA O 00 o 
fs.r^fOfOr*-r>-o»-» c o 
CM rH X CM X CM X rH 


*-h m cm x m rH m n 

rH rH X X 


\XOOVmOOCMrHrH 

rOOvCMvoCMCMr^r^ 

vo^OXCMCMCMCMCM 


A A A A A A 

mCJHHvtH 

mm rH 


O' CM X X N ov 
rH vo X X X 00 
m m vo vo x rH 


vaaaaaaaaaaaa 

<mvOXvOC 7 WOvor^CMCMCMrH 

rH CM rH CM rH rH 


/A A A A A A A 

HK 1 CM vo CM CM CM 

HHNH 


moCMr^cMCMmmfMmooovooorH cm o rH m m 

K)(M(SJf / )>OK)ps.NvOO\fOOff'N^ >0 CM VO N 

CO X CM X CMCM CM CM m M) vo rH X OO CO CM X CM CM CM 
>4 rH rH rH 


! ^X HfOHH K) CM r» ’ 

^SSSSESSSSSSoSSSSSS^SN^kNvrrtWO^vOvOrHCNJMCNIC'lNC'JinvOVO'O^r-NO^NINJCNJCJrN 

. . ..... I i — < • — I rH i — 4 rH rH H rH rH 


AN AN AN AS AN AN AN AN AN AN AN AN 


CM CM CM 

vo vO vO 


o o o o 
oooo 


mx m vo 

rH H H H 
tOKIfOK) 


oooo o 
coco oo oo co 


vo vo VO vo vo vo 


o 

o o 

oooo 

o 

oo o 

o 

o o 

oooo 

o 

ooo 


cm m 

XtOvON 

00 

coo o 

rH 


rH rH rH rH 

rH 

O O rH 

m 

mm 

m m m m 

m 

m m m 

a 

o o 

oooo 

o 

OOO 

CM 

CM CM 

CM CM CM CM 

CM 

VO vo vo 

m 

m m 

m m m m 

m 

mm m 

so 

vO vo 

vo vO vo vo 

VO 

sO vO vO 


oooooooo oo oo oo 

oooooooo oo oo oo 

rH cm m x m vo r- oo co ov oh cm m 

rH rH rH rH rH rH rH rH OO rH rH HH 

K)K)K)K)K)fOKiro mm m m mm 


vO vO vO vO vO vO vO vO OO OO OO 

mmmmmmmm xx xx xx 

vO vO vo vO vO vo vO *0 vo vO vO vO vO VO 


ot ov m 

UJ m rH 

CO rH 


cm O' o> x o 

o noiooo m 

m rH rH 


no co x co r^oocMx vo cm oo o 

so vo cm cm vo m co m co 

rH rH rH CM rH 


r^xovoor^mxN x oo cm ov com 
xcofN. mcM vo vo m x h-x 

CM rH rH *“» r-i r+ 


<r ITinDN 

CO 

o 

rH 

CM 

m x 

m vo f-N oo 

in m m in 

in 

Ov 

OV 

Ov 

Ov Ov 

O' O' O O' 

CM CM CM CM 

CM 

CM 

CM 

CM 

CM CM 

CM CM CM CM 


cm m x m vo r>. «o ov ov o rncM mx 

mmmmmmmm vo r- £[>- £-£- 

mmmmmmmm mm mm mm 


93 



/A 

A A 

y-A 

A 

A 

^A *r> 

Pa 

CM CM 

to 

to 

o 

00 CM 

m 


CM 


CM 

sT 

Ay* 

Ay* AT 

Ay* 


Ay* 

Ay* <k> 

VO 

NT rH 

to 

H 

SO 

O' o 

O 

os in 

SO 

sO 

rH 

sr m 

00 

sT to 

pA 

CM 

PA 

0 O rH 


iH rH 


rH 


rH 


CMsT 

mc\i 


o to 
NT CTs 

oo -H 


CM 


0 "T sr ms 


ST rH S3- K> 
h- (M O' O' 
COOvJ-H 


O' 

m 

r>- 


/-> a» 
NT 00 


v£> LO 
tO vO 
CM <\l 


*~s tr\ ^ 

r-» to m CM 

rH lO SO 


IO 

r- 

so 


s 3 - so m 
Pa o o 
co oo so 


r> ^ 
rH CM 
to 


O' rH 
m CM 

o 


o 

CM 


ST 

p- 

co 


VO O 
r-H IT) 
fA- H 


SO 

to 

CM 


(OtOvJs 

mrHtOrH 


0 O' io o 
ST ST sO m 
00 00 tO rH 


AS AS /—\ 

r> (Mo i—i pa 

lO Nj- 


a: 

UJ 

« 


m 

o o 
so CO 


O CM 

in cm 

—I <M 


Pa N i— < I'. 

H HNr-< 


o 

to 

SO 


to so m 

SO *— < rH 

pa h- m 


/A ArA«A 

*H 


y-A 

' 3 - rH CM 0 O 

O' 

m rH 

to CM 

'T 

Ay* Ay S-T Ay* 

Ay* 

to SO 

Ay* W 

m 

Ay* A> 

OOrH sT sT 

to 

sT so 

0 O rH 

NT Os O' Cts 

SO 

PA O 

sT CM 

SO so ST sT 

Pa 

00 oo 

SO O 


^ ^ A, AAt AVAN AS 

•nn to sr to h*. cm 0 

• CM so CM to i-h CM 


0000 rH tO 00 O' rH <X 

oooo to rssr nno 

PaPa 00 so so CM O sT 


c 

o 

o 

a; 

nJ 

CQ 


m ST rH pA CM to 

sr 


cm O' O' O' m O' 

so m vt n sr rs 
iONOOfsJoOH 


i—i O' 
to rH 

oo m 


AAA 

A A AS As 

AS 

A A A A A 

o to O 

rH O' m rH 

to 

ST HH rH CM 

•h m »h 

rH tO 

rH 

rH in rH ST rH 

SrfS/W 

Ayr A> Ay* A>* 

Ay* 

Ay* Ay Ay> Ay* Ay* 

to CM CM 

m oo m m 

PA. 

SO tO sO CM O' 

pA SO pA 

rH ST SO SO 

OO 

O Pa rH Pa >r 

SO in sT 

m so to to 

m 

OO so fA. \T 00 


co 

CM 


O 

m 


AAAAAAAAAA 

cMtomocMtomoocMto 

toto in to CM CM 


o tOtOA-tsOvTcOCMsOsom 
CM so NO to O IV. sr K> A-« so 

r-v fv. fv 00 00 00 so ST CM fv. to 


pa- 

pa 

rv 


a 

Q 

13 

o 

IS 

a 

& 

o 

a; 

O 


JS 

1= 

a; 

a? 

** 

O 

ns 

co 

nj 

OJ 

GO 


oc 

UJ 


A AA A A AA ^ rA 

A AS 

AS y-A A 

sT sr oo so rH in m 
m sr rH 

CM so 

m 00 rH 

to 

A/ S-* A* A* W W A* 

A> Ay* 

A_y S_> Ay* 

m m CO so rH ST pA 

SO O' 

CM CM rH 

rH O NT tO O' ST tO 

pA. CM 

SO Pa fA. 

m -so so cm so m sr 

ST sT 

in sT CM 


A AAA A A A 

in 0 "j- to m cm o 

tO w — I rH 


rH 00 CM CM CM CM sj 

pa m cm cm o cm sr 
cm m to to oo cm in 


AAAAAAAAAAA A AS 
OOIOHOOOOH\T\XHO o 
AHtO tOHMfsjH rH ^ 

W'XWWSyiWWWu'S^ W Sty w 

0 rH o in in oo oo ct' cMn rH cm is. 

h to to o o cm sx pv in 'o O' o m 

in 00 so M> SO CM so CM Pa to >>o ©0 CM 


A AAAAAAAA AAAAA 

' 3 , sTNsTsOs 3 -HinkOMoOrvK )0 

SO aH tO to rH rH 


S 3 - N N O tO O' O' tO rH Os rH O' O' 00 

s 3 -CO 0 OtO|s.NNO'O'NO'tf|iOm 

mimn'oso'T'cr»HsocM'orv|v.tn 


AAA 

00 to to 
m sr 


**S AS As 
tO Pa H 


a A A A 

O' to CM O 
CM 


rH OO 00 

r-m m 

cm m m 


H\JO\ 
O' sj” OO 
SO so — l 


ST fv. CM in 
OS OO to fA. 
to m CM CM 


^ ^ ^ ^ ^ ^ AA aaaaaaaaaaaaaaa AA _ _ ■ 

saas!8s~-~»-"—ssasssaa-3s-2 — ssSSsssFS™ 

S58SSlip3|S3^^^i2^3iS3SRxSi§ii!g3S^S 

^ rH rH rH rH 


QJ 

o 

o 

o o o o o o 

DOODOO 

O O O 
O O O 

o o o 
o o o 

O O 

o o 

ooooooooooooo 

OOCMOcOsOsTCMOcOsOsTCM 

O 

o 

oooooooooooo 

COSOCMOOOSO'O-CMOOOSO'T 

o 

CM 

oooo 
O 0 O so ST 

yO 

ns 

Eh 

3: 

<*o 

o 

T 

sr m so pa oo oo 

rH rH rH «H rH O 

to to to to to to 

O' O rH 
O rH rH 

to to to 

CM to ST 

rH rH rH 

to to to 

315 

316 

P>-COpAOOcOOsOHCMCMtOsTin 

HHOOOOHHHHHHrt 

tototototototototototototo 

316 . 

'OrApAOOcOO'OrHCMCMtOsT 

rHHOOOOHHHHrtrt 

tOtOtOtOtOtOtOtOtOtOtOtO 

m 

rH 

to 

so so pA. 00 

rH rH rH rH 
to to to to 


0 

1 

o o o o o o 

OO O O O ST 

o o o 

sT ST sT 

oo o 

'T sT sT 

o o 
sr sr 

ooooooooooooo 

STSTOOOOOOOOOOOOOOOOOOOOOO 

o 

00 

oooooooooooo 

COCOHHHHHHHHHH 

o 

rH 

OOOO 

HHHH 


LAT 

vT sr ST sT sT sT 
sO so SO so sO so 

ST ST ST 
SO SO so 

ST ST ST 
so so sO 

ST ST 
SO so 

sTSTvTsTSTSTSTsTsTSTsTSTsT 
SO so so so SO sO so so so sO SO so sO 

ST 

SO 

'T'Tmmmmmmmmmm 

SO SO sO SO SO sO S) so so sO so sQ 

m 

so 

m m mm 

s£> so so so 


to 

t- 

Q_ 






UJ 

to 

2 

M3 rH CM so Pa O 
SO tO rH rH iH Pa 
CM rH 

to to Os 

rH CTS 

O O to 

CM rH 00 
rH CM rH 

189 

120 

tOMnrHmtOtOmMOHNrH 

rH rH »0 ST O' *0 tO H m 

CM CM rH 


r-v cm »-h to o rH sr o so oo to so ia. 

Ca H CM to rH rH tO "O Pa O' SO 

CM tO 


sa- 

in 

CM 


O' S3- m CM 

to CM 


QC 

UJ 

to 

2 E 

=> 


2 

PH 


in so pa oo c?s os 

pAfA. PA pA pA O 

to to to to to sT 


o rH cm to sr m 


sT ST NX NT sT \T 


so |a. 

H r— l 

sT ST 


00 O' o rH CM to *T m SO pA CO Os o rH 

rHrHininminininininioinso >o 

sT'TsrsT'TsTsTsTsT'TsT'TsT sT 


CM tO O rH CM tO »T m so Pa oo O' O 
'O'OOOOOOOOOOO rH 

'T'Tininininintninminin in 


rH CM tO sT 


in in min 


94 



CM rH SO 
CM 


O' LO so 

sT VO PO 


As A A a A 
vT\3-OCON 
CM PO sT CM PO 


in Nino sx 
HHvTK) CVJ 


/s rs ^ 

CM \X CM so 

sX PO rH CM 


C?s sX m rH sX 

\T rs VO M ^ 
ooeoi^ovr 


osHino sx 

sT PO so m C- 

00 00 PO -H 00 


Os 00 O' PO 
vTOON O' 

00 r**» CM H 


/*> a rs ^s rs r* a a> ^ 

H0O SOHOsOMC r- P"* 

rH LO rH PO PO PO rH 


so so H O' O' sT m VO 

cd o po sx c- h- so »-h 

oo to oo oo cm oo po r** 


a nrs Arvn^sn rs ^ 
PO 00 CM SO CM CM >3“ 00 PO M) 

K> CM sX CM PO >3- rH sx 

w SW W W W S-/ SW W W SW 

SO «0 pH O' pH sT pH sJ" SO O' 

o oo os c- cm O' os p*- rH m 

oo r*. so cm o st so oo n r- 


so m sr so 

J H N 


CO rH so O 

sX O' pO in 

so so CM rH 


A A AS AS AS AS AS A AS A AS A AS 

CMoOOf-iintjCM'OsTCM«-tin ©0 

rH rH m K) 'O fO H H r-» 

vy sy sy s> vy v»y »*> s> vy s> w 

00 «H«K)vOHOsvOOIOff'N 

eocofO'TNfOO'rs.HOSHi^rs 

rs.fs.oososocMsoCMr-'^mr-p- 


As yA AS A AAA An A A A A 

hO po so CMsTsommCMOsTsT «H 

CM rH rH sX M> CM m 

ww w w sw sw s> w w w w W sy 

o-po in oo oo so rH m po >o so cm so 

rH SO o sX SX PO PO SO P"- I""- rH ps. P- 

\X fs. so so so CM 00 PO M> so P". *sT m 


AS A A A A A AAA A A A 
sOOfsHvO LO sO sO O m PO rH 

CM rHsJ- CM CM f — rH J) 

VO m sT OO PO 00 sT rH rH vX PO pO 

© O sT sX Os oO O' fO PO ps ps fs. 

00 so m SO rH ps- sT CO 0O 00 so sO 


AAAAAAAAAAAAAAAA As AS a AAA, AA a AAAA 

rH OS CO SO s3“ Sj CM CM OO CM SO sT so P^ so o P- rH sOCOPOCMsOOSOOOPOsTsT PO 

K) CM pn rH CM PO P^- PO CM CM POsTPOsTPOPOrHPO rH 

s>wsy^s>wws>wsysys/ws>syw ws> sysy w w w w WWW ww ^ 

m cm cm m m o *> po po so cm cm o m oo sxos vtoons^hosoopoon in 

op^sosooohOsrosr^-por^r^cTsrHin © *-t 'Xoococooscsr>.po<>r^co a* 

in cm in*m sosososorH'ocMsT’X'Tmin cm m mp-mmrH'OCM'osxsocM m 

rH r-4 rH rH rH 


sT sT P- P- sT m sX 00 CM sO 00 CD sT m OO 

rH CM >-* CMrHPOm <T 


os CM SO m rH P- sX o rH in r^. o rH OS so 

h ro o o o m s3 m cp so co »o po cm r- 

sr o oo so m cm m rH <«o po m so oo sx m 


ors.r^.s 3 -s 3 T^r^OOr-'OsPOOsoOOOO'Os(MO'inr^r^' 3 ‘OOC 30 'OOOSO'COr^rHOOOin' 3 'CM<MrHr^OOrHs 3 -r^O'OOOOsO' 

SKEnS5««SSnSosnj?MCMcoooorHr^inin'3*'osoinpOK)osCMCMCMp-.rs.CMin»oososooin'osoCMc7soor-^<>pooocM 
~ S S S S " S ^ ^ S 2 >r CJ (vl ” " co US w PJ CVI to r. r! r. VO CVI ^ vr vr CVJ CVJ « o (SJ « as ^ (VJ r, r, o ST ui ^ vr SO CVJ vr 

^ rH rH HrH rH rH rH 


ooooooooo 


o o o o o o o o o o o o o o 

CMOOO'sOsS’OCMO'O'sTCMOOO SO 


OOOOOOOOOsOOOOOOO 

' 3 ‘CMOOO'O'O’CMOinCMOCOsosTCMO 


P's. OO CO O' O rH CM CM PO 
OOOOHHHHH 
POPOPOPOPOPOPOPOPO 


sj- m SO 'O P* 00 o P- oo O' o rH CM CM PO 

rH HHHHHNOOOHHHH rH 
ro K> PO PO PO PO PO PO PO PO PO PO PO PO PO 


srm'O'Ops-ooo'OrHrs.ooooo'OHCsj 

HHHHHHHNNOOOOHHH 

POPOPOPOPOPOPOPOPOPOPOPOPOPOPOPO 


o o © o o o 
OO so ST CM O OO 


cm po sr m so so 

rH rH rH rH rH rH 

PO PO PO PO po PO 


ooooooooo 

sX'T'T'XsX'T'T'XsX 


oooooooooo o o o 

sXsXsXsXsXsxp-p.p-p-p-p-p* 


oooooooooooooooo 

i^f^p.p.p.p.p.p.r^ooo©ooo 


ininminminminm 

sO'D'O'O'O'O'O'O'O 


m m in m in m m m m m m m m m m 

SO sO so so sO sO sO sOsOsO s O'OsO'0 so 


inininininmininm so-so so so so so so 

'OsOsO'O'OsOsO'O'OsO'OsOsO'OsOsO 


DODOO 

o o o o o 


sO SO sO sO SO 
sO sO sO sO sO 


PTS 



u 

sJ-'J'sflOCPH'OlOsO 

m 

UJ 

cm rH m in m so sr 

PS. 

OQ 

s: 

ZD 

Z 

rH PO 

sX 


in cm cm so m oo po CM rs- so o po m rH po o po m o p- sx po oo cm ih p*. o m p- po cm 

Is. POOSfHPOPO rH HNinOP' rH U"| P- CM sX SO H CM rH CM PO O rH 


SO fs- sj- sO sT 

CM m so sX Os PO 

CM rH rH 


23 o rH CM PO s3" in so rs. 00 Os O rH CM PO vTvOOHPO'T m so p» oo 

z m m m m m m m m m m sosososososoooooooo o 

in m m m m m in m m m m m m m m m so so so so so so so so 


OSOrHCMPOs 3 ‘in\OOOCDrHCMPO' 3 *msO 

OrHrHrHrHrHiHrHrHintnminminin 

SO'O'^SOSO'OSOSOSO'O'O'OSOSOSOSO 


95 



c 

o 

O 

0) 

Cfl 

cC 

C3 


c$ 

Q 

3 

o 

s 

a 

03 

&Q 

O 

0> 

O 

c 

JS 

c 

a> 

a; 

J-. 

o 



PH 

PH 


PH 

PH 

PH 


PH PH 


vr 

rH 


NO 

nT 

NO 


f-v pH 







rH 


pH 


hp 

HP 


HP 

HP 

HP 


V/ HP 


o 

CM 


ON 

o 

o 


NO (M 


CM 

CM 


rH 

uo 

uo 


pH CM 


CO 

vr 


nT 

pH 

pH 


r-v >T 



pH 


rH 

pH 

pH 


i—4 

• 


PH 

PH 


PH 

PN 

PH 

PH 

A PH 


m 

pH 


rH 

nT 

r-v 

pH 

NO NO 








pH 



HP 

HP 


HP 

HP 

HP 

HP 

HP HP 


NO 

ST 


NO 

rH 

P-4 

NT 

CM On 


pH 

NO 


f-v 

(VI 

CM 

rn 

no r-v 


c- 

pH 


NO 

O 

o 

00 

NO NO 



pH 


pH 

rH 

pH 


pH 


PH PH 

PH 


PH 


PH 

PH PH PH 

AAA 

PH 

rv in 

UO 


nT 

NO 

rH 

m m i-v 

NO rH NO 

00 

(VI CM 

pH 




CM 

NO rH NO 

CM 

H 









Ql 

NP HP 

HP 


HP 

HP 

HP 

HP HP HP 

HP VP HP 


NO On 

60 


NO 

pH 

pH 

rH pH ST 

CM On O 

QC 

pH pH 

uo 


NO 

ON 

ON 

no no r-v 

r-v r-v uo 

LU 

fH.NO 

UO 


NO 

NO 

NO 

60 00 00 

NT NO CM 

A 




pH 




pH pH 

X. 









=> 









z 









VP 

PH PH PH 

PH 

PH 

PH 

PH PNPNPS 

PH P\|H APlPHiHAA 


CM NO NO rH 

NT 

NO 

no r-. co cm 

l-v ON NO O O CM NT CM nT 

> 

UJ 

NO 



pH 

rH ST pH H 

CM NT pH nT CM NO rH 

a: 

HP HP HP 


HP 

sp 

HP HP HP 


HP s^r VP HP HP HP 

HP HP w 


NO UO ON in 

NO 

nT 

ON rH O' nT 

eo r-v r-v no no no oo no on 


fv. pH UO rH 

NO 

O 

nT O f-v nT 

oo r-v f-v r-v r -- ph p-i no uo 


m m r-v m 

r-v 

CM 

oo m cm m 

r-v fn r-v so no r-v nt no r-v 





pH 

rH 



H 


-PH PH ^N 

PH 

PH 

PH 

rs/-srv 

PH PS PH PH PH PN PH PH PH PH 



NO rH O' O' 

CM 

CM 

rHNOrHeOO>OOf-vONOr-v(MNOOOinUO 


CM NO 






rH NO >T NT NO pH pH 

H 


HP HP HP 


HP 

HP 

NPVPNP 

V> HP 

HP HP HP HP HP HP SP HP 

WWW 


CM CM UO 60 

ON 

CM 

WNNOHONOONSOOvCCONCOOlTi 


r-v |w rH rH 

pH 

NO 

\3-infOONcocooofOfOrvHHMifiH 


ntntuo vr 

uo 

O 

NOCMCMinoMf-vininNONDinNONOCMCMin 


M5 NO 
NO 


sT so 
r-. no 
« o no 


>r 

o 

(VI 


CO 

60 

rv 


sTCMf-v 

*-h vr 


vr cor-, 
vr eo i-'. 
icr^N 


/N^N^^N/N ^N 

60 no nt as no r-v 


O'OHHI's h- 

l-v loo O r-v eo 

(si r-i m in r-v uo 


o 

m 

CM 


O pH cm 


no no 
r^r- 

NO NO 


P> P> PH PH AS PA 

CM CM O' 00 vT pH 

m 


pH nT cm cm sr m 
NO P-- 'O NO NO "O 

00 60 no NO pH 


PH PH PH PH PH AS PH 

N^Ohh CM nTCM 

CM (VI 


OHO^ no CM no 
CM NO rH CM HfMH 
«0 60 NO nT r- NT C- 


PH PH AS PH AS pH, ^H AS AS PH AS PH PH PH PH AS pH PH AS PH PH pn PH PH PH 

►Olf\M , HNO'N'O l OHO'M , HHO'M , OvTvOO'(iOeO'0«nH 

rH NO CM pH UO rH NO NO nT NO 


NO 

Ch 

m 


PH PH PH PH PS PH P 

rH NO (VI NO rH NO C 


p- P> PH PH PH PH PH PH PH PH PH PH PS PH PH 

oinCMrH.\oeoNOoc3NovrCMeor-ipH 

CM rH NO NO NO UO CM rH 

HPSAWSPHPSPHPHPSPWHPSPHPN^HP 

M’NOOM‘'TNOOO'00'CoO(MM 

ONOOON'sTvreoNONOCMrHNoCMrHrH. 

(ViONX^UOUOlTl'ONO'O'MrNTCMNONT 


CM 

o 

60 


1 PH PH PH PH PH PS PH PH PH PH PH PH PH PH PH PH PH PH PH PH P 

| pH-ONO|H.OOUOpHNOsONOCMNOCMlOrP>rHCNJp-^N3-UOr 

I CM NO \T no r-l NO CM NO r-l 


OOOOHpHHlfUMMv£(>OO^NV)M-oaNONOCvONM-M'(OCNf\JNN'OvOfOO'rHNN<>CN'OvOON<ONTCT'M->OcOOa'n-»rH 

CM(ViCHr^or^NONOON3-NONoo'mNOvroNooooiH.(vicMrHrHNOrHiriO\3*soOrHNo^<MOirioocOfHfH. r -« r HrH.rHr^«>CMiOiOCNo 

CMCMCMCMNoCMuoin6oeorHp^^(MCMuo'«TCMCMLnsrvTNTrHrHvrrH.ooinin60NON060NOCMCMCM<viuoir\NOvrr-<NOrHCM(v»rH.rH.rHNo 


</} 

S3 

a> 

C/2 

CD 

a 

3 

as 

H 


z 

QC 

O 

u 


OOOOONOOOOOOO 

nOnTCMOMOnTCMOPOnonTCM 


ooooooooo 

eONOVrCMOeONONTCSJ 


OOnOn^OOOOOOOOOOOOOOO 

ONN)HM-CslO«0^'3 > O00'0MTMO00'ONr 


o 

CM 


O On O' O' O' O' 
OrHinriOfs. 


MS 0 

r-eooNOHNOMOoooNOH 

pHpHpHCJCMOOOOOHpH 

NONONONONONONONONONONONO 

CM 

pH 

NO 

CJNONTinNONOr-vCOO 

pH pH pH pH pH pH pH pH pH 

NONONONONONONONONO 

OOCMNOvor-vOOCOONOCMCMNOsTinvOsOI-vOO 

CMCMCMCMOOOOOpHpHpHpHpHpHpHpHpHpH 

NONONONONONONONONONONONONONONONONONONO 

On 

pH 

NO 

O O pH NO no nT 
CM CM CM CM CM CM 
NO NO NO NO NO NO 

z 

0 

—l 

1 

oooooooooooo 

O O O O O NO NO NO NO NO NO NO 

O 

NO 

OOOOOOOOO 

NONONONONONONONONO 

ooooooooooooooooooo 

NONONONOsOSOnOnOVOVOnO'O'O'O'OvOnOnOnO 

o 

-o 

O O O O O O 
SO NO SO NO NO NO 


NO NO NO NO nO no so nO NO NO no no 
NO NO *0 NO NO 'O <0 NO NO NO NO NO 


NO 

NO 


nOnOnOnonOnOnOnonO 
nO nO'OnOnO'O'OnOnO 


nOnOnO'OMJ'O'OnOndnonOnOnOnOnononOnOno 

sOnOnOnOnOnOnOnOnONOnOnOnonOnOnOnonOnO 


NO 

NO 


NO NO M3 NO NO NO 
NO NO NO NO NO NO 


C/0 

h- 

a. 

Of 

tu 

« 


60NoON60CMinpN.voONr-«iNOio 

m sT O CM rH CM CM ON 


MThHON (VIOlHiO 
rH.rH.NorHONNorvrH.-H 
CM rH pH pH 


ohocviNONOKOOrMNTyoKjooNONONomcovrrH. 
NO rH HH NO rH m Pi rn CM CM O rH to 


pH NT CMrH H ITTr 

in Nj- pH pH 


QC 

LU 


CQ 

j: 

=> 


NOnJ-UOnOCOOnOhCMNOnTIO 
nO NO NO NO nOOnOOOOOO 
NONONONONONorHrH.|H.r,|H.rH. 


NO fv- 60 On o pH (VJ NO sy in 

O O O O pH rH pH pH pH pH 

fS. fH. f-v |H. |H. fH. fv. JH. fn. IW 


NorvO'OONOHCJN>M - NONeooNOH(viro\f 

pHpHrHCMvrinininmcnmmmmNONONONONo 

N |H- fN. fN f-v f-v fN- fN. fH |H. fv f-v fs. f-v f-v |-v f-v ^ fv. 


uo 

NO 

r-v 


nOChCOOpHCM 
NO NO NO f- IH. f-v 

f-v fH. fv. rv JH. f-s. 


96 


590( 




<rv 


rv 


00 00 vT rH oo 

rH m 

00 

VO 

to CM 

CM CM 

rH 

CM 

to 







v> v/ v> v> vy 

vy \/ 

w 

vy 

vy vy 

vr m v3" m O' 

\T O' 

O 

rH 

vr m 

'O'O'O'ON 

O' r^- 

CM 

to 

vO vO 

HvJr-tVTtO 

vT CM 

OO 

oo 

rH vr 

hhhhh 

f — i rH 



i — < «-h 





G 

O 

O 


in 

»- 

o_ 

a 

UJ 

CQ 

2 : 

=> 

z 

> 

UJ 

a: 


as 

PO 

as 

aS 

Q 

13 

y 

IS 

a 

& 

o 

<u 

O 

c 

<D 

CD 

Sm 

O 



r\ 

yv yv 

yvA yv 


to oo n- to oo 

n- rH 

vO 

CM 

CM O' 

VO rH 

r- O' to 

vr 


rH rH CM 

vy vy vy vy vy 

v/ vy 

V/ 

vy vy 

vy vy vy 


mcM m cm to 

VO MD 

r- 

00 

CM 

rH CM M) 


to CM to CM O' 

to to 

r- 

OO 

vO 

tO CM rH 


O vr o vr CM 

CM CM 
H H 

n- 


vO 

oo vr r- 

rH 

y*N 

rH rH rH rH rH 
rsyv Ayv a 

y-s y~v 

yv yv 

y^ rv yv 

AnnyvAA 

vr 

O' O' vO «0 O' 

fv CM rH 

r-o 

CM rH 

rv o vr 

CM CM H 

m VO 00 CM rH O' 
VT M3 rH tO 

vy 

CM rH 

vy v/vyvy vy 

<^y> W W 

w v> 

vy vy vy 

wvy vy v> vy vy 

O 

0O O' CM O' CD 

to to S3- 

r- vo 

O' O O 

CM O' O' vr vT vo 

CM 

oo r- no r- m 

O' O' O 

r- n- 

rH CM tO 

^ Nfvrv.0 h 

CO 

to so to cm 

*— » i— t CM 

HHH 

r-m 

VO OO vo 

vo to to oo m r^- 

rH rH 


H rH rH 


r^- vr o vr 
CNJ 


CD O CM rH 
Ifl^KIvO 
rlHOH 


to c\j o r-* m «o 


tOQ'NNfON 

to HH 


hhoonco^ 

oooooocoeo 

iniAN^INlA 


inKlNHOtO 
vOHHl/1 (SI 


CD K) VT tO CM CM 
too'r^r^h.N 
vo CM 00 vO <vT <t 


^ ^ rv /"v ^ ^ 


m 

r- 


r.^ 

00O'»n>J->ON in m CM CM vO CM m CM 

sT cm cm m r-t m »-* 


in 

CD 

M 3 


«OhhhN 
vT CM CM O m 
vO © O vO O 


Ov Ov <0 vO Ov 00 

(OWNNHH 

cm cm mm vo vr 


COOIOHOONO 

r-tinr^MDvovovrov 

vTCMsOMTCMCMmin 


(M 

o 

oo 


CM 

CM 

to 


r— I CM 00 CM in 

rH rH to 


O' o' cm cm r- 

mtoKiN 

CM VO o CD 


to rH r*- 

moosiniAto 

VO rH 

rH 

rH rH 

CM tO CM 

CM tO 

to 

vy vy vy 

vy v> vy vy vy vy 

vy vy 

vy 

vo O' O 

ooNsorvvorv 

1—4 rH 

vT 

r- rH CM 

0O 'O to oo to CD 

to to 

n- 

in VO oo 

rv. so to m to cm 

rH rH rH 

OO OO 

oo 


^ ^*s rv <— s rA y— v ^ ^v ^ iTV av 

rHin'jfor^O'OCMmvr\T^r-( 

to CM to rH CM CM 'J" 


HH(Mj'Noo'OooN r )inmoo to 

OOHHvOHOHOtOvOtOfO f"- 

ininvovovO'TCMvrcMOMrsovo md 


to CM ■O’ 
to 


rH vO 0"0 


^ y— V y-V /—>. y 

vovy com c 
cm cm m m cm M> * 


v^vrv^/v/ 


s y— v ^v y-\ /— s rv rv av rv 


IO tO M3 

r — •— • 

m m vo 


VO O' 
r-H m 
VO M3 


O' o 

m o' 

vO vT 


ONh-vJ 

ovin m vr 
vj cm cm m 


SS^S^SSSSSS^SSwinintoin^'TM-'OCMCM^incMCMCMinin^cM^cMinvTvcr'T'OCM 


as 

CD 

m 

dD 

CD 

IS 

173 

H 


03 

UJ 

z 

Qd 

O 

o 


to 

o 

z 

0 

—I 

1 

h- 

< 


OOOOOOOOOOOCDOO 

NTMOW'O'fCMOCO'O'S’CMOCO 


OOOOVOVO'O'O'O'OOOOOOOOOOOOOOOOOO 


'on.coooO'OHCMCMtova-invo'O 

OOOOOrHHHHHHHHH 

to»otototototo^ototototototo 


rH hhnnncm(sjn(Moooo( 


OOOOOOOOOOOOCDO 

0 'OvOvO'C 7 ' 0 'OvO' 0 'C 7 vO'OvO'< 7 v 


S S O' Ov <> CM CM CMCMCMCM^CMCMCMCMCMCMCMCM CM CM CM CM CM 


VO'O'O'O'O'OVOVOVO'O'O'O'O'O 

vO vO vO vO vO MD vO vO vO vO vO vO *0 *0 


MDM3vcvovovovovovoi^r^r^f^r>-r^r^f^c^r^^c^^ r ^c^r^h~r>- 

VO'OVOVOVOVOVOVOVOVOVOVO'OVO'O'O'O'OVO'O'OVOVO'OVOVOVO 


o 

o 

CD O O O 

o 

m 

o un om 

r- 

r- 

00 OO O' O' 

i— 4 

rH 

rH rH rH rH 

to 

to 

to to to to 

o 

o 

o o o o 

CM 

CM 

CM CM CM CM 


r^- 

r^r^r-r'* 

VO 

vO 

VO VO VO so 


<o 

t- 

0 - 

a: 

UJ 

CO 


MjCMvotommcovor^CMCMCMoor- 
to VT CM to to CM to to OO m Vf rH 


VO 

to 




O' 

O' 


CM COST CM O' 
VO sOh.H'J" 


oc 

UJ 

pa 

2 : 

=5 

z 


0 'OHCMtO'ffmM)N «0 0 'OH(NJ 
OvCOOOOOOOOOHHH 

^coeoooeotooooooooofiooococo 


^invor-oocTv^CMto^rinvor-ooovtovrinvor-wovo^rMto^ 

_,—» r -.—i^-iH(MCMCMmiriinmminvo' 0 'Ovovovovon*-rv-r*-r s *-n : 

Ooooooooeooooooooooooooooocoooooooooooooeooooooooooooo 


m 

r^ 

oo 


VO 

r^. 

00 


r*. co O' 

r*- r- tv. 

oo co oo 



97 


088 


c 

o 

O 


^3 

DQ 


cs 

Q 

3 

y 

IS 

§• 

& 

o 

cd 

O 

c 

3 

<d 

CD 

s* 

O 


CO 

I— 

CL 


or 

LU 


s: 

2 


> 

tu 

DC 


v3* 

e*. 

oo 


no 

h*. 

VO 


f^CM 


O' CN 

r-- r- 


VO vo 

rH rH 

r^- r-- 


vo 

eg 


eg 

eg 

NO 


V/ vy 

eg vo 
eg no 
no eg 


Ov O 

e- m 
eg pi 


<^n r> rv 

o vr eg 


ST rH 


eg 

oj eg eg 

o o NO 


/~v /■>, rv r“\ 

eg no eg vj- oo 'O eg 


ov r*. eg ej 

eg vt r- r^- 

v3- lo <r sr 


h eeovf >o 

eg sr sr no r-s. no 

eg oo oo uo r>- eg 


r>. /s 

LO ‘sO 


eg eg 

NO NO 
o o 


in vT LO 

eg 


IflNN 
vo pn. 



rs rv 

y> rv 




m so 

O O' 




eg eg 

►o eg 




w vy 

vy w 




sr ej 

ST ej 




so eg 

so eg 




rH vT 

rH ST 




rH H 

rH rH 




/**s 

rr, rs 




rH|v. 

me 




no eg 

eg eg 




w vy 

vy w 




me 

m O' 




no e* 

NO 




o NO 

O NO 




rH rH 

rH rH 





■NS r-v 


«r> 


in no r-* 

00 O' 


hv «gj 


eg no eg 

eg 


NO rH 


vy sy w 

w vy 


w w 


eg eg vo 

00 vo 


rH VT 


vO SO NO 

00 NO 


no e- 


VO SO NO 

f'. NO 


CO CO 


rH 

rH 



r-V 


yv rv 

<r> 

^rvrvrv 

Ov 

sT CO Nv. 

o o> 

rH 

VOHPsH 


ih(M 

vt eg 

rH 

eg rH 

w 

sy vy vy 

vy w 

vy 

wwwvy 

O 

rH *0* NO 

g-io 

*o 

O rH NOrH 

eg 

O sT O 

VT O' 

CO 

NO NO pv O' 

OO 

m m eg 

meg 

r- 

VO 00 VO h- 


rH 

rH 



/•v «— v 

Arv/v 

y% rs 

rirv^v 

<r> y> r> r*v i 

NO NO NO O vo 

m o 

emmemH'croi 

eg eg <—i nt eg 

K> 

no eg no 

NO rH 


o> 

r- 

eg 


vo 

H 

e- 


o 

m 


CT> rH CO 00 O CO o 

HOHH1T rH m 


m 

NO rH 

rH 


w 

vy vy 

vT 

NO O 

NO 

sr m 

CO 

eg rH 




a 

r Ji 

7i 

CD 

m 

C£> 

JD 

3 

c3 


0C 

O 

o 

2 

00 

o 

2 

0 

—I 

1 

h- 

< 


ooooooooooooo 

oLnomoinoininmomo 


oooooooooooooo 

ininoinuioinoinoinomo 


oorH^-tegegNONOinvoe-r^co 

egegegegegegegegegoooo 

NONONONONONONONONO»nNONOfO 


ooooooooooooo 

egegegegegegegegeggrvTvjvT 


ooovooHegegNONOvrva-uoinvo 

OOHHHHHHHriHHHH 

nonononononokjnonononononono 


oooooooooooooo 

v 3 ’vT’a'v 3 -v 3 -s 3 ->T N a'v 3 '«d'vtf t Nj*>J-> 3 - 


r^f^r^i^r^r^r^e-.r^r^r-r-e- 

s O' 0 ^^g}\ONCvO' 0 «OvOvOvO 


r^r^r^f^e-r^r-e-h-e-e-e-e-e-. 

v©v©'0'g>vOvOvO\©vO'OvOvOv©v© 


oooooooooooooo OOOOOO 
oinoinoinmoinoinoinmoinoinoin 


g’vfteieg) vo pvN*o*oovovoHHNcgioK)g , iesosfv.ooco 

r? ri| rj ^ ^ r? hhhhhhcnj(\J(NJ(\I(\j(\joj(s|(\iooooo 


r>- 

vo 


nononononononononononononononononononono 


vS-va*\|-<r^grv 3 'srgrvTva->TvTv 3 -sTvOvOvOvOvo 

^' , ** s *e s *e , v , e ,i *^ s *f ,,, »r ,ta r s -ewr s -e-e , ^r ,,, ^r^e«.e^e«.e^rN. 

vovovovovovovosovovovovovovovovovovovovo 


CO 

I— 

0l 

DC 

LU 

CQ 


eo vj- vr o rH eg no eg ►o sj- ^o va* eg o rH no st rH sa- eg e- eg rH \j- eg rH co 

Nocg no no rH eg eg m rH eg vt eg no no so vo vo 

NO 


eg 

NO 

NO 


^vovof^oovTNominNoegNoegegegrHrHNoavvo 
in in no rH «o no no no c\igg eg _« 


CQ 

z: 


•HrgNO\Tuovor-.coegvj-invoh- 

0O 00 00 00 CO 00 00 OO O' NO NO NO NO 

oooocoooooooooooeoovovovov 


oo o rH eg >r m vo oo ov o rH eg no vt 

no vr gr va- sj* vj- vj- xj- ^ nj- m in m m m 

Ov OV O' Ov Ov C?v Ov Ov Ov 0% Ov Ov <Jv Ov Ov 


mvON«OvO(\JK)vnnvONC)OOHvTlflveNOO 

iniemioiegj'O'ovovovovo^NNHHHHH 

Ov O' O' Ov O' O' CN Ov O' O' O' ON Ov O Ov O O O O O 


98 


530( 6) 77<K 10) 


Table 6. Seasat Greenland Geographical Data Base (Cont.) 

BIN NUMBER NUMBER PTS LAT-LONG SW CORNER REV (NUMBER PTS) 
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Table 6. Seasat Greenland Geographical Data Base (Cont.) 

BIN NUMBER NUMBER PTS LAT-LONG SW CORNER REV (NUMBER PTS) 

1537 3 68.80 328.00 228 ( 3 ) 

1538 1 68.80 328.50 228 C 1 ) 

1540 1 68.80 329.50 515 ( 1 ) 
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Table 6. Seasat Greenland Geographical Data Base (Cont.) 

BIN NUMBER NUMBER PTS LAT-LONG SW CORNER REV (NUMBER PTS) 

1751 15 69.40 515.00 243 ( 15 ) 
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Table 6. Seasat Greenland Geographical Data Base (Cont.) 

BIN NUMBER NUMBER PTS LAT-LONG SW CORNER REV (NUMBER PTS) 

2149 63 69.90 319.20 286 ( 26 ) 530 1 20 ) 774 ( 17 ) 
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Table 6. Seasat Greenland Geographical Data Base (Cont.) 

BIN NUMBER NUMBER PTS LAT-LONG SW CORNER REV (NUMBER PTS) 

2477 6 70.20 330.39 788( 6) 
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Table 6. Seasat Greenland Geographical Data Base (Cont.) 

BIN NUMBER NUMBER PTS IAT-L0NG SW CORNER REV (NUMBER PTS) 

2626 19 70. AO 310.00 A32( 11) 1264( 8) 
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Table 6. Seasat Greenland Geographical Data Base (Cont.) 

BIN NUMBER NUMBER PTS LAT-LONG SW CORNER REV (NUMBER PTS) 

3053 72 70.80 320.79 901( 18) 645( 23) 819 C 8) 896 C 23) 


Os 

CO 

oo 


POCsl 

rH CM 


K) 

CM 


W V-C 


SO OO 
ST -H 

oo o 


oo 


AN /^l AN ^ 

sT PO PO \T CM no 
CM CM H r — 4 r -4 


Os O Os t© CO nT 
rH 0 O OO *-H 0 O PO 
OO CO OO O SO O 


AS 

AN 

AN 

SO 


PO 


rH 

CM 

w 

w 

s_r 

rH 

to 

O 

ON 

t*N 

SO 

oo 

oo 

h* 


AS AN AN AS AN AN AN 

AS AS 

AN 

AN 

CM so O' rH lA PO sT 

POO 

PO 

PO 

rH HH rH 

«H 

CM 

CM 

Nysysywsrs/sr 

w w 


o 

00 0 O LOO O O O 

NN 

CM 

NO 

sr st rs. cm on cm on 

rH rH 

PO 


00 CO OOO so o so 

r^r>* 

CA 



O' 

CO 

oo 


SO 

«T 

00 


i-Hr^STPOK)sTPOr-l 

NrHHMNMNH 


HHflOTHHfs. 

CM rH 


sT O' oo r- r^h-tM 


rH PO PO 
rH CM CM 


mosooeosTr^sTsT 

NTrHOOOOsTrHPOPO 

vOflOsO'O'I'f oo 


PO 

ia 


N^s/wsrs/UN> 

w w 

s/srsrs/sr 

w 

s«* s> 

CM CM 0 O rH so o NO 

ON OS 

O so SO O' SO 

NO 

SOLA 

PO PO NT On stf- On NT 

OO oo 

^OHNlflN 

sT 

NT sT 

00 00 00 OO sT no nT 

nT sT 

r^ir»fN.ioi^ 

00 

00 so 


CM'T'TsTCsJrHsTPO IT) O sT 

H CM i —4 CM H r -4 »—4 CM CM 


rHrHCMsOlAPONTrHsXPO rH h* IA CM PO K> M> 00 PO h- sT PO ia CM po r*- 
H r-t CM CM CM hhHNhN CM HtSlH CM «— * 


oo oo oo nt h- o r^ 

H H H H O H 
SO M> so <f <r CM sT 


ro 

ST 

CM 


O si sj 
nONN 

>t a-a- 


s-c w v> s> s> w s> w v> s> 

inifUVJNNMflOOO 
OOPOstTTNMfOfO 
00 OO OO so SO SO 00 sT sT \T 


K>K)O"O*0CM(MininP0ososcjsu*unr^. 

SNCOHCOPOK)|s-NOHHH\JTH 

sTsTsrmcMimnioinoooooooososo'T 


^s ^S AN, AN A ' AS an A /A rs A As A A AS AS AS As AS AS A A A A A AN AS AN AN A An an AN As A AN AS AN AN AN AN AN AN AN AN an AN AN AN AN AN 

KOO'Mmp r )fosr> / )HsrNOK)soirioNr-4fsjHrjHHoorNrN«-<o«HC\nooNrcsiNa , K)ino^K)^sTorHK)N)OHiriiri 

CM rH CM CM CM CM CM ( rN rH aH rH rH rH CM rH rH <M CM rH rH CM rH CM CM CM rH CM rH 


HHvIOOK)K)K)oO' 0 ' 00 "f K)OCy|00«0OK)Ne0HH00HHHHO^'O'C»0C>0NCbW(M«0«0O'^UV^«000C0HK) 
OOvTOONNN'0'OCOflOlDrsrsCPK)rHHPOOHinNOsOlOKlK)K)K , )POCOOO«rsCJ(M(MNNOO«imOrs.OHHHOrN 
sJsITNCMhhh sTsrCMfNJCMPslOfNsOfs.N'O'O'OCMlfHOirisOsO'O'O'lHHHCtMMCSICMCMCNiNlOlOineO'O sOs^slH 


0s0s0'0s0'0‘0'0'0'0'0'0s0'0s0'0'0s0s0‘0'0s000000000000000000000000000'0'0s0' 

A 4 ir» 0 'K)r^A-imOsPOrNrHmONPOrNONK)inCMN-ir»NTCMOOO'OsTOOCM'X»Os 3 -oOCM' 00000 (MsoO'T« 0 (MsoOs 3 -fN.A-»iOO' 


HHH(MCMK)POK)sJ> 3 -knUHnvO'ONOrONlvCN'a , ONOOHfJNrO»ONl<sf tflOlOvcsO'ONMOlOCOO'^OOOHHH 
NNCMCSJC\J(MN(MCM(\i(M(MCMNNNK)Nl>OK»»OOOiHHHHHHHHHHHHHHHHHHHiHHrtN(M(MNMC\J 
POPOPOPOrOPOPOPOPOPOPOPOPOPOPOPOPOPOPOPOPOPOPOtOPOPOPOPOPOPOrOPOPOPOPOPOPOPOPOPOPOrOrOPOPOPOPOPOPOPOPO 


OOOOOOOOOOOOOOOOOQO o oooooooooooooooocaooooocjooo o OO P o o 

OOoOOOOOOOOOOOcOOOOOOOOOOOOOOOOOOOOOOOOOOOCTsOsOsOnO'O'O'O'OsOsOnOsOsOiOOsOnO'O'O'OsO'OsOsO'OsO'O'OsO'O' 

ooooooooooooooooooooooooooooooooooooocooooooooooaooo 


MOrv 
o ON so 


sOO>JiHfs.HsJPOON'OtOOC\|HNA(CMHHeO««OPOOONiHCMCOOO(MMOsOO»nPOCSJNrK)cONO 
so r^. fN. NO »H CM CM CM PO CM rH rH rH rH CM CM PO sT CO no sT rH IA so CM sT Os so IA sT nT CM rH <T 


OO ST |H CM CO 

po cm co r^~ o 


TiriNorNfioooHNfO 

LOIOIOIOIOIO'O'OSOSO 

oooooooooo 

POPOPOPOPOPOPOPOPOPO 


sTIOsor^OO 

so SO SO SO so 

a a ooo 

PO PO PO PO PO 


rH ION© 

a- a. oo oo 
o o o o 

PO PO PO PO 


POinCMsTsOOOOCMPOsTlOvor^eOOsO 

CTNOrHCMCMCMPOPOPOPOPOPOPOPOPONT 

OOHHHHHHrlHHHHHHH 

POPOPOPOPOPOPOPOPOPOPOPOPOPOPOPO 


rHCMPONTirtsorNOOON 

sTnTsTnTn 3 ’nTnTsT , <T 

HHHHHHHHH 

fOPOPOPOPOPOPOPOPO 


OHNPOvf UtNON 

ialapaialaialaia 


POPOPOPOPOPOPOPO 


118 


322.39 1 7 3 C 23) 917 C 29) 695( 3) 889( 20) 1018( 17) 1039( 21) 
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Table 6. Seasat Greenland Geographical Data Base (Cont.) 

BIN NUMBER NUMBER PTS LAT-LONG SW CORNER REV (NUMBER PTS) 

387 A 7 71.60 329.19 645C 7) 
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Table 6. Seasat Greenland Geographical Data Base (Cont.) 

BIN NUMBER NUMBER PTS LAT-LONG $W CORNER REV (NUMBER PTS) 

4252 514 72.00 320.40 158( 21) 20H 21) 244C 21) 287C 21) 445C 21) 488C 21) 

6 Q3( 21) 617 ( 21) 646( 21) 689( 21) 775( 21) 804( 21) 

847 ( 21) 890( 21) 1019C 20) 
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Table 7. Seasat Geo-referenced Data Base Header Description 

FILE 1: GEO-REFERENCED DATA BASE HEADER RECORD 

Record Format: One logical record corresponds to one physical record 

Blocksize: 480 Bytes 


Bvtes 

FORTRAN 

Variable 

Type 

Description 


1*4 

Number of latitude rows in the data base (56) 

5-8 

1*4 

Northwestern-most latitude of data base in degrees North (x 10 5 ) 
(7210000) 

9-12 

1*4 

Northwestern-most longitude of data base in degrees East (x 10 5 ) 
(30000000) 

13-16 

1*4 

Southeastern-most latitude of data base in degrees North (x 10 5 ) 
(5990000) 

17-20 

1*4 

Southeastern-most longitude of data base in degrees East (x 10 5 ) 
(34000000) 

21-244 

1*4 

Width of each latitude row in degrees (x 10 5 ), starting with the 
southernmost row. This is dimensioned by the number of latitude 
rows in the data base. 

245-468 

1*4 

The number of longitude divisions in each latitude row, starting 
with the southernmost row. This is dimensioned by the number 
of latitude rows in the data base. 

469-472 

1*4 

Logical record in data base at which directory starts. 

473-476 

1*4 

Size of the data base, including the directory, in blocks. 

477-480 

1*4 

Status word for altimetry data. 


Bits 

value 

Description 

0-23 

0 

Unused 

24 

1 

Slope correction applied 


0 

Slope correction not applied 

25 

1 

Orbit adjustment applied 


0 

Orbit adjustment not applied 

26 

1 

Solid tides removed 


0 

Solid tides not removed 
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Table 7. Seasat Geo-referenced Data Base Header Description (Cont.) 


(477-480 Cont.) 


Bits 

Value 

Descriotion 

27 

1 

Re tracking correction applied 


0 

Retracking correction not applied 

28 

1 

Center of gravity bias applied 


0 

Center of gravity bias not applied 

29 

1 

Tropospheric correction applied 


0 

Tropospheric correction not applied 

30 

1 

Ionospheric correction applied 


0 

Ionospheric correction not applied 

31 

1 

Time bias applied 


0 

Time bias not applied 
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Table 8. Seasat Geo-referenced Data Base Description 


FILE 2: GEO-REFERENCED DATA BASE 

Record Format: 595 logical records correspond to one physical record 

Blocksize: 19040 Bytes 


Subgroup 1: 

One logical record for each bin containing data 

Bvtes 

FORTRAN 

Variable 

Type 

Description 

1-4 

1*4 

Indicates the number of logical records which follow which are 
located in the bin 

5-32 


Unused 

Subgroup 2: 

One logical record for each data point in the bin 

Bvtes 

FORTRAN 

Variable 

Type 

Description 

1-4 

1*4 

North latitude of datum point in degrees (x 10 ) 

5-8 

1*4 

East longitude of datum point in degrees (x 10 ) 

9-12 

1*4 

Surface height relative to the ellipsoid in cm. 

13-16 

1*4 

Height sigma, arbitrary value of 1.0 m used (x 10 ) 

17-18 

1*2 

Rev number 

19-20 

1*2 

Used for temporary flags when gridding the data 

21-24 

1*4 

Orbit adjustment in meters (x 10 5 ) (-999999999 if unavailable) 

25-28 

1*4 

RMS of orbit adjustment in meters (x 10 5 ) (-999999999 if 
unavailable) 

29-32 

1*4 

Slope correction in meters (x 10 5 ) (-999999999 if unavailable) 


NOTE: Subgroups 1 and 2 are repeated for as many bins with data. 
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Table 8. Seasat Geo-referenced Data Base Description (Cont.) 


Subgroup 3: 

Directory 


Bvtes 

FORTRAN 

Variable 

Type 

Description 

1 — 4 

1*4 

Record number 

5-8 

1*4 

Record number 

9-12 

1*4 

Record number 

13-16 

1*4 

Record number 

17-20 

1*4 

Record number 

21-24 

1*4 

Record number 

25-28 

1*4 

Record number 

29-32 

1*4 

Record number 

NOTE: 

The directory contains as many 3 
the record locations of all bins. 


at which data for bin 1 starts 
at which data for bin 2 starts 
at which data for bin 3 starts 
at which data for bin 4 starts 
at which data for bin 5 starts 
at which data for bin 6 starts 
at which data for bin 7 starts 
at which data for bin 8 starts 

2-byte logical records as necessary to designate 


Table 9. Elevation Grid Header Description 


FILE 4: 

ELEVATION GRID HEADER RECORD 

Record Format: One logical record corresponds to one physical record 

Blocksize: 80 Bytes 

Bvtes 

FORTRAN 

Variable 

Type 

Description 

1-4 

1*4 

Number of latitude increments in the grid for a non-polar 
stereographic grid (140) 

5-8 

1*4 

Number of longitude increments in the grid for a non-polar 
stereographic grid (152) 

9-12 

1*4 

Starting north latitude of grid in degrees North (x 10 6 ) (this will 
be approximate for a polar stereographic grid) (50000000) 

13-16 

1*4 

Starting east longitude of grid in degrees East (x 10 6 ) (this will 
be approximate for a polar stereographic grid) (300000000) 

17-20 

1*4 

Ending north latitude of grid in degrees North (x 10 6 ) (this will 
be approximate for a polar stereographic grid) (73000000) 

21-24 

1*4 

Ending east longitude of grid in degrees East (x 10 6 ) (this will be 
approximate for a polar stereographic grid) (340000000) 

25-28 

1*4 

Status word for data used to generate grid. A zero in any bit 
position indicates that the correction is not applied. 


0 


31 

Bits 

Value 

Description 

0-23 


Unused 

24 

1 

Slope correction applied 


0 

Slope correction not applied 

2b 

1 

Orbit adjustment applied 


0 

Orbit adjustment not applied 

26 

1 

Solid tides removed 


0 

Solid tides not removed 

27 

1 

Re tracking correction applied 


0 

Retracking correction not applied 

28 

1 

Center of gravity bias applied 


0 

Center of gravity bias not applied 

29 

1 

Tropospheric correction applied 


0 

Tropospheric correction not applied 

30 

1 

Ionospheric correction applied 


0 

Ionospheric correction not applied 

31 

1 

Time bias applied 


0 

Time bias not applied 
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Table 9. Elevation Grid Header Description (Cont.) 


Bvtes 

FORTRAN 

Variable 

Type 

Description 

29-32 

1*4 

Polar stereographic grid size conversion and scaling factor from 
half-inch grids on projection plane to the desired grid size (x 10 ) 
(1650000) 

33-36 

1*4 

The number of grids of desired size from the pole to the equator 
based on the grid size conversion and scaling factor (x 10 ) 
(608754894) 

37—40 

1*4 

Latitude of the map perimeter in degrees North (x 1 0 6 ) 
(500000000) 

41-14 

1*4 

Greenwich orientation in degrees (x 10 6 ) (450000000) 

45-48 

1*4 

Polar stereographic switch (1) 

=0, grid has constant increment in latitude and longitude 
= 1, grid is in polar stereographic projection 

49-52 

1*4 

Number of I-axis divisions to the extent of the map perimeter 
(445) 

53-56 

1*4 

Number of J-axis divisions to the extent of the map perimeter 
(445) 

57-60 

1*4 

J coordinate of the projected pole (223) 

61-64 

1*4 

I coordinate of the projected pole (223) 

65-68 

1*4 

Minimum J index of the grid (166) 

69-72 

1*4 

Maximum J index of the grid (317) 

73-76 

1*4 

Minimum I index of the grid (305) 

77-80 

1*4 

Maximum I index of the grid (444) 
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FILE 

Bytes 

1-4 

5-8 

9-12 

13-16 

17-20 

21-24 

25-28 

29-52 

53-76 

77-80 

81-84 

85-88 

89-92 

93-96 

97-180 

NOTE: 


5: 


Table 10. Elevation Grid Description 
ELEVATION GRID DATA RECORD 

Record Format: 10 logical records correspond to one physical record 

Blocksize: 1800 Bytes 

FORTRAN 

Variable 

—Type Description 


Condition number of the matrix used in the least-squares solution 
to the function (x 10°) 

Capsize in degrees latitude - radius from grid location defining 
area from which data was used to define grid (x 10 6 ) 

North latitude of grid point in degrees (x 10 6 ) 

East longitude of grid point in degrees (x 10 6 ) 

Height values of the grid at location relative to sea level in meters 
(x 10°) 

Number of data values that were used to calculate grid value 
Number of parameters used to define function, NPT, (equals 0, 


Six gridding function coefficients. If NPT is < 6 then the rest of 
the coefficients are initialized to zero, (x 10 5 ) 

Set of null coefficients associated with any negligible singular 
values (see SVD reference). If NPT is < 6 then rest of coefficients 
are initialized to zero (x 10°) 

Distance in km from grid locations to closest data point (x 10 6 ) 

Nortl^ latitude of closest data point to grid location in degrees 
(x 10) 


** 4 East longitude of closest data point to grid location in degrees 

(x 1 0°) & 

** 4 Height associated with closest data point to grid location in 

meters (x 10°) 

1*4 Standard deviation of the data with respect to the gridding 

function in meters (x 10°) & 6 

'* 4 Correlation matrix from solution. This is a symmetrical 6X6 

matrix so only the upper triangular portion is stored. The order 
of storage is elements 1-6 are the first row elements, 7-11 
columns 2-6 of second row etc. (x 10 5 ) 

Ten of the above-mentioned 180-byle logical records make up one block of data. 
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16. Abstract 

The data processing methods and ice data products derived from Seasat radar 
altimeter measurements over the Greenland ice sheet and surrounding sea ice 
are documented in this first volume of a series. The corrections derived 
and applied to the Seasat radar altimeter data over ice are described in 
detail, including the editing and retracking algorithm to correct for 
height errors caused by lags in the automatic range tracking circuit. The 
methods for radial adjustment of the orbits and estimation of the slope-induced 
errors are given. The various levels of ice data sets are described in this 
report, but the user is referred to Volumes 2 (Greenland) and 4 (Antarctica) 
for more detailed descriptions of the gridded elevation data sets and the 
geo-referenced data bases. 
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